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"A ce compte, toutes les sciences ne seraient que des applications inconscientes 
du calcul des probabilites ; condamner ce calcul, ce serait condamner la science 
tout entiere." 

"From this point of view all the sciences would only be unconscious applications 
of the calculus of probabilities. And if this calculus be condemned, then the 
whole of the sciences must also be condemned." 
Henri Poincare, 1914, La Science et I'hypothese 

Abstract Solar gravity modes (or g modes) - oscillations of the solar interior for which 
buoyancy acts as the restoring force - have the potential to provide unprecedented 
inference on the structure and dynamics of the solar core, inference that is not possible 
with the well observed acoustic modes (or p modes). The high amplitude of the g-mode 
eigenfunctions in the core and the evanesence of the modes in the convection zone make 
the modes particularly sensitive to the physical and dynamical conditions in the core. 
Owing to the existence of the convection zone, the g modes have very low amplitudes 
at photospheric levels, which makes the modes extremely hard to detect. In this paper, 
we review the current state of play regarding attempts to detect g modes. We review 
the theory of g modes, including theoretical estimation of the g-mode frequencies, 
amplitudes and damping rates. Then we go on to discuss the techniques that have 
been used to try to detect g modes. We review results in the literature, and finish by 
looking to the future, and the potential advances that can be made - from both data 
and data-analysis perspectives - to give unambiguous detections of individual g modes. 
The review ends by concluding that, at the time of writing, there is indeed a consensus 
amongst the authors that there is currently no undisputed detection of solar g modes. 

Keywords Sun ■ theory ■ data analysis • g modes 



1 Introduction 



The detection of solar oscillations bv lLeighton et alj (jl962h marked the start of a mem- 



orable era for solar physics, and the birth of a new field, helioseismology. The oscillatory 
motions detecte d in the solar atmosphere had typical periods of about five minutes, and 
iDeubnerl (|l975t ) subsequently verified observationally that the oscillations were the vis- 



ible manifest ation of acoust i c (p) modes trapped in cavities in the solar interior. Later 



in the 1970s, Iciaverie et al.l ( 19791 ) showed that the Sun support ed truly global modes 



of oscillation. The low-angular-degree (low-i) modes observed bv lClaverie et all (1979) 



penetrated the solar cor e, and had therefo re opened a window on the structure of the 
deep interior of the Sun. iGrec et"aL I l|l98d) subsequently made long-duration observa- 



tions of the Sun from the South Pole, providing the very first clean power spectrum 
free from day-night interruptions. 

While inference on the internal structure of the Sun made using high-quality data on 
the p modes has gone from strength to strength over the intervening 30 years, attempts 
to detect g modes have proven altogether more elusive. This is in stark contrast with 
the ubiquitous detection of stellar g modes acro ss the Hertzsprung-Rus sell diagram 



such as in white dwarfs ( Winget fc KeplerlboOct ). in subdwarfs B stars ( Green et al.1 
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2003), in B stars l|WaelkenJl99lh and in 7 Dor stars (jAerts et al.lll99Sh . The detection, 



and subsequent exploitation, of data on the solar g modes would provide much more 
precise inference on the structure and dynamics of the solar core than is possible with 
p modes. This is because the g modes are trapped in cavities in the radiative interior, 
and they may have much higher displacement amplitude in the radiative zone / core 
than their acoustic counterparts. However, this high sensitivity to the core properties 
comes at a price. Because the g modes are trapped beneath the convective envelope 
they have very low amplitudes at the photospheric level where the perturbations due 
to the oscillations are detected; hence their elusive nature. 

The first attempt to detect long-period oscillations date back to the mid 1970s, 
e.g.. ISevernvi et"al] 1 19761) . Further cl aims, from ground-b a sed D oppler velocity obser- 



vations, were made in the 1980s, e.g. JPelache fc Scherrerl l|l983h . It was realized very 



early on in the development of helioseismology as a field that a space mission would 
be ideally suited not only for observing the p modes but also for detecting g modes 
(in principle providing the long-term stability needed to detect the low-frequency g 
modes). The DISCCp mission was proposed by iBonnet et al.l ( 198ll ) primarily as an 



irradiance mission to which helioseismology was added later on. The initial design was 
subsequently modified to accommodate severa l in-situ instruments and spectrographs, 
and finally matured into the SOHQj mission (jPomingo et al.ll 19951 ). 



One of the goals of the helioseismology instruments on SOHO was to detect g 
modes. Following its launch in 1995 December, initial analyses of the helioseismol- 
ogy data failed to uncover evidence for the modes. By 1997, it was realized that the 
searches would benefit greatly from an internationally coordinated effort comprised of 
the principal observational and theoretical scientists involved in g-mode research. This 
led to the formation of the Phoebus Group, which has been working since to use all 
available high-quality helioseismic data, and to develop the analysis techniques that 
are required, to detect g modes. 

In this review we collect and summarize the knowledge accumulated by the Phoe- 
bus group over more than a decade of g-mode research. Our review is broken down as 
follows. In Section 2 we review the theoretical properties of the solar oscillations. [A 
complete prese ntation of the prop e rties of solar and stellar oscillations may be found 



complete presentation 01 trie properties ot solar and stellar oscillations may be round 
for example in Unno et al.l d 19891 ) , IChristensen-Dalsgaard fc Berthomieul () 199 ll ) , and 



IChristensen-Dalsgaard ( 2002al ).] We also consider the accuracy of the theoretical com- 



putations of the g mode frequencies, and detail the various factors that contribute to 
the uncertainties (i.e., the accuracy of the oscillation computations themselves, and the 
sensitivity of the frequencies to the structure and dynamics of the solar models) . 

In Section 3 we discuss the excitation and damping of the g modes, again from a 
theoretical perspective. We consider the different physical mechanisms that can poten- 
tially contribute to the stability of g modes, and present theoretical estimates of the 
amplitudes of g modes excited stochastically by convection in the convective envelope. 

In Section 4, we turn to the observations. We begin by considering the different de- 
tection techniques that are available to observers, and discuss in detail the constituent 
elements required to "build" a front-to-end analysis algorithm. Our presentation quite 
naturally includes an in-depth discussion of frequency-domain analysis, and hypothe- 
sis testing (frequentist versus Bayesian). We finish this section by reviewing (from the 
literature) attempts to apply these techniques to detect g modes. 

1 Dual Spectral Irradiance and Solar Constant Orbiter 

2 SOlar and Hcliospheric Observatory 
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Finally, in Section 5, we conclude with some remarks on the current state of play. 
Our view is that, at the time of writing, there is no unambiguous evidence for the 
detection of solar g modes, although some of the most recent theoretical predictions 
suggest a positive result may be attainable in the near future. We therefore finish by 
looking to the future prospects for finally detecting g modes. 

2 Eigenfrequencies for the solar model 

The Sun can show oscillations whose restoring forces are due either to compressibility or 
to buoyancy. Compressibility gives rise to acoustic waves: resonant pressure (p) modes 
in the high-frequency range (typically about 3000 fiRz). The effects of buoyancy can 
give rise to internal gravity (g) modes, which have very low frequencies, and are the 
focus of this review. The so-called f modes (surface gravity waves) and mixed modes 
(having dual p- and g-mode characteristics) occur at frequencies around the transition 
between g modes and p modes. 

2.1 Properties of gravity modes 
2.1.1 General properties (adiabatic) 

The Sun, as a mass of compressible self-gravitating gas, can oscillate around its equi- 
librium state. The small-amplitude oscillations displayed by the Sun may be described 
by spherical harmonic functions, Y [ m (9,(j)), of degree I and order m where (9, (f>) are 
spherical polar co-ordinates, 6 being colatitude (i.e., the angle from the polar axis) 
and 4> longitude. For example, for a singlet, the perturbation of a thermodynamical 
quantity like the Eulerian pressure perturbation p' can be written: 

p' =p(r)Yt m (e,4>)e- 2 ™^ mt , (1) 

where v n ^ m is the cyclic oscillation frequency (which may be related to the angular 
frequency by to — 2ttv). The displacement vector of the material from its equilibrium 
position £ may be written, in the frame of polar spherical coordinates, as: 

I = ^(r)Yr,L-^ h (r)^Yr,imL-%(r)^j e -*™^\ (2) 

£ r and fih being related to the radial and horizontal components of the displacement and 
L — y/lQ, + 1). Each mode is characterized by three numbers, n, I, m: the radial order n 
is essentially the number of nodefl of the radial displacement eigenfunction from centre 
to surface in the radial direction, the degree I is related to the horizontal wavelength 
of a mode, and corresponds to the number of nodal lines on the solar surface, and the 
magnitude of the azimuthal o rder m, where m satisfies — I < m < I, is twice the num- 
ber of nodes in lon gitude (e.g. lUnno et al.lll989l : Ichristensen-Dalsgaard fc Berthomieul 
ll99ll : lGoughlll993h . If the Sun were spherically symmetric (no rotation and no mag- 
netic fields) the frequency v n i m of a mode would not depend on the azimuthal order 
m. 

3 At the limit of low and high frequency, the radial order n is essentially the number of nodes 
of the radial displacement in a solar model. But the classification of modes, i.e. the assignment 
of a radial order to them, may be mor e complicated depending on the frequency range (case of 
mixed modes) and on the model (e.g. lEckartlll960i IScuflairejll974l : IOsakilll975l : lTakatall2UO§) . 
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Fig. 1 Propagation diagram for a reference solar model. The solid and dashed curves represent, 
respectively, the real values u;+/2-7r and u)~ /2n corresponding to zero values of K 2 (r), and are 
plotted as a function of the radius, for 1=1. The hatched areas indicate p-mode and g-modc 
propagation regions for this degree. The full horizontal lines indicate the propagation regions 
of three modes of degree 1=1: a p mode of frequency v ~ 600/iHz (p3), a g mode of frequency 
v ~ 150^tHz (g3), and a mixed mode of frequency v ~ 280/iHz (pi). The corresponding 
cigcnfunctions of these modes are shown in Figure [2] 



The solar oscillations have two restoring forces, compressibility and buoyancy, char- 
acterized respectively by the sound speed c and the buoyancy (Brunt- Vaisala) frequency 
TV: 



2 

C = 



N*=g 

P 



1 dp 
7ip dr 



1 dp 
p dr 



(3) 



Here, g, p and p are, respectively, the gravity acceleration, the pressure, and the density 
of the solar model, r is the radius, and 71 is the adiabatic exponent. In order to have 
an insight into the relation between the oscillations and the structure of the model, one 
may consider a simplified wave equation , which provides an approximate description 
the oscillations (|Deubner fc Gougblll984h . If the wavelength is much less that the solar 
radius, the local effects of spherical geometry on the dynamics can be ignored. This 
equation is then derived under the local wave-like approximation JWKF4J, and under 
the Cowling approximation, i.e. where the perturbation to the gravitational potential 
by the waves is neglected. The equation is 



d 2 <P 
dr 2 



+ K(r)& = 



(4) 



with 



K(r) 



s? 



(5) 



The quantity & is related to the displacement of the wave, via & = cp" 
uj c and H are, respectively, the Lamb and the acoustic cut-off angular frequencies, and 



4 JW KB stands for Jeffreys- Wentzel-Kramers-Brillouin {Jeffreys I925|). See also iGouehl 
in the asteroseismic context. 
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Fig. 2 Kinetic energy density (i.e pr 2 |£| 2 normalized by the inertia X n as function of the 
radius for m odes g3, g2, gi, pi, P2 and p3 of degree 1=1, as computed for a reference solar 
model. From lProvost et all l|2000l V 



the density scale height: 

In this plane-wave approximation, k r represents a local radial wave number, which may 
be associated with a local horizontal wave number, = L/r. In the wave equation, the 
term in iV 2 is the signature of the buoyancy, which acts as restoring force only in the 
radiative regions (where N 2 > 0) . The stratification in density close to the solar surface 
induces an acoustic cut-off angular frequency, ui c , which corresponds to an upper limit 
to the frequencies of trapped, radial modes (1 = 0). 

Equation @ shows that a mode will oscillate as a function of r in regions where 
K{f) > 0. Such regions are called regions of propagation. In contrast, a mode will be 
evanescent, i.e. its energy will vary exponentially, in regions where K(r) < 0. K(r) is 
positive either for ui > (uj + ,uj~) or w < (ui + ,w~) , the angular frequencies uj + and 
uj~ corresponding to the zeros of K(r) . cu + and lu~ are close to Si and N respec- 
tively, except in the external layers where u> + is close to u> c and lu~ is very small. The 
frequencies uj + and u>~ are plotted, as functions of the radius, in the so-called propaga- 
tion diagram, for a solar model (Figure [TJ. The horizontal lines represent modes; they 
are continuous in the propagating regions and dotted in the evanescent regions. There 
are two possible regions of propagation, according the frequency of the modes: a high 
frequency region (u> > Si,ui c ) corresponding to pressure modes (or p modes) resulting 
from the compressibility, and a low frequency region (u) < N) corresponding to internal 
gravity modes caused by the buoyancy. In between, some modes, called mixed modes, 
can experience a substantial restoring force from both the gradient of pressure and 
buoyancy. 

In the high-frequency range, which corresponds to observed solar oscillations, the 
modes are acoustic modes, and their properties depend mainly on the variation of the 
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sound speed c(r). Their radial wave number fc r is 




Thus these modes obey the classical dispersion relation fc 2 + fcj 2 ~ and the fre- 
quencies increase with the radial order n. Owing to the increase of the sound speed 
towards the solar interior, the p modes are refracted at the level where lo = S[. Thus, 
they are trapped in a cavity between that level and the surface. The lower the degree 
of the mode, the more deeply it penetrates the solar interior. An example of a trapping 
region is indicated in Figure Q] for the p mode with n = 3 and 1 = 1. Figure [5] shows 
how the amplitudes of various modes vary within the solar interior. 

In the low frequency range, we are dealing with gravity modes (g modes), for 
which buoyancy in radiative zones (iV 2 positive) provides the restoring force. Their 
radial wave number is given by: 

N 

k T ~k h —. (8) 

The properties of these modes therefore depend principally on N, and the frequencies 
decrease with increasing radial order n. Because g modes oscillate in regions where 
their frequencies are smaller than the Brunt- Vaisala frequency, g modes are confined 
in the solar radiative zone and core (see Figure [T] and Figure [2}. Low-degree gravity 
modes have much larger displacements in the horizontal than the vertical, the opposite 
of what is shown by p modes. 

The g modes are attenuated in the convection zone wit h a factor proportional to r^ 7 , 
r C z being the radius at the base of the convection zone (|Christensen- Dalsgaardiri980l ; 
iBerthomieu fc Provost]|l99d) . so that low-degree g modes may be easier to detect than 
their higher-degree cousins. We shall see in the next subsection that the frequency 
spectrum shown by the g modes is manifested as a pattern of roughly equidistant 
periods for modes of given degree I. 

In an intermediate frequency range between 200 and 400 ^iHz, there are low radial- 
order g modes, f modes (surface gravity oscillations that have no radial nodes) and pi 
modes (a mode with one radial node) . Some of the modes, called mixed modes, have a 
pronounced mixed character, i.e. they have amplitude both in the central layers, like g 
modes, and in the external layers, like p modes (e.g. mode pi (I = 1) in Figure [2l ). The 
properties of mixed modes have been studied in detail by IProvost et al.l (|2000l ) . The 
mixed character of these modes appears very clearly when we look at their inertiae, 
I n i . The inertia may be written 

ln,i = [ Pt-tdV, (9) 
Jv 

where the eigenfunctions £ are normalized at the surface of the Sun and where the 
integral is taken over the whole volume of the Sun. Figure [3] plots the normalised 
kinetic energy of the modes, which is proportional to the mode inertia, as a function 
of frequency. Under the assumption of equipartition of the energy in the modes, the 
surface amplitudes woul d be approximately inversely proportional to the square root 
of the mode energy (e.g. IBerthomieu fc Provostll 19901) . 

For g modes with frequencies less than 200 /^Hz, the lower the degree, the smaller is 
the energy, and the higher is the surface amplitude. Around 280 /^Hz, modes of mixed 
character have smaller energies and therefore higher surface amplitudes than modes 
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200 400 600 800 

Fig. 3 Logarithm of the normalized energy ui I n ; of low- frequency g and mixed modes 
(o pen symbols) and lo w-frequency p modes (filled symbols) for a reference solar model Ml 
of lProvost et al, | feOOCft , plotted as a function of frequency, for modes of degree Z=0, 1, 2, 3, 4, 
5 and 6. (The normalization is taken assuming that the modes have equal amplitudes at the 
photospheric level, where the temperature equals the effective temperature level). Note the 
transition from p to g modes around 450 /iHz, and the existence of a set of modes of mixed 
charact er around 280 f iHz, having lower energies than modes in the neighbouring frequency 
regions Provost et al. (2000). 




Fig. 4 Relative contribution of the buoyancy restoring force to the frequency, /C2/(/Ci +AC2). 
fC± and K.2 are defin ed in Eq. dlOI . for the same modes, model, and symbols in Figurc[3] From 
iProvost etHTl 12000ft . 
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adjacent in frequency. We discuss the excitation and damping mechanisms of g modes 
in detail in Section [3] 

The fr equencies of the os cillation modes may be expressed using a variational prin- 
ciple (e.g. lUnno et alj Il989h . Neglecting surface terms and the perturbation of the 
gravitational potential, we have: 

(2™) 2 ~ J- f ( dr + J- f N^pr'dr S (10) 

where the Eulerian pressure perturbation, p, is given by p = pru L The quantities 
K,\ and K.2 correspond, respectively, to the contributions to the frequency of pressure 
and buoyancy forces. 

Figure [4] shows the variation of K,ij{K,\ + IC2) as a function of the frequency, for 
the modes of a reference solar model. As expected, for p modes there is a dominant 
contribution to the frequency from the "acoustic" term K\, while for g modes it is the 
the "buoyancy" term, IC2, which dominates. Mixed modes, between 200 and 450 /iHz, 
have frequencies defined by significant contributions from both terms. 



2.1.2 Asymptotic approximation 

Low-degree g modes with c yclic frequencies les s than i/„ i < 200/^Hz may b e described 
by an asymptotic relation (|Vandakurovlll96Sl ; iTassoull [l98ol ; lOlverl Il956h . For small 



frequency, the second-order asymptotic expression for the period P n ; (P n .i = l/^ra,/) 
can be represented by: 

P n ^P n , l = P ^(n + l/2-\ + ,) + ^^^, (11) 



with P a = — = , (12) 



Here, ■& is a ph ase factor sensitive to the properties o f the layers lying below the 
convection zone |Provost fc Berthomieu1ll986l : [Eliij|l986l '). If A^ 2 is assumed to vanish 



dr 



2 



proportionally to [r C z — r) v for some power p at the base of the c onvection zone, •& tends 
at low - frequency to a constant — 0.5/(p + 2) (see also for details iBerthomieu fc Provostl 
Il99lh . For the standard solar model, a linear behaviour of A^ 2 may be assumed and •& 
tends to —1/6. The second-order coefficient V\ depends on the B runt- Vaisala, frequency. 



V2 d epends in a more complicated way on the stratification ( Berthomieu fc Provostl 

S). 

In sum, the above tells us that the frequencies of g modes are related closely to 
the Brunt- Vaisala frequency, particularly through Pq and Vi . Typical values of Pq , Vi, 
and V2 for a reference solar model are, Pq ~ 35 to 36min, V\ ~ 0.4 and V2 ~ 5.7. 

The periods of g modes of a given degree I are proportional, in the first order, 
to Pq/L. The separation in period between modes adjacent in n is therefore almost 
equidistant. In the low-frequency asymptotic limit, the relative frequency difference of 
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two modes of the same degree, I, and same radial order, n, from two different models 
is equal to: 

$v n ,l sp nd 5P 



(14) 



Deta i led comparisons of the numerical periods and the asymptotic periods (|Provost fc Berthomieul 
1 19861 ; IChristensen-Dalsgaard fc Berthomieul Il99ll ) have allowed checks to be made of 
the validity of the asymptotic formula. A least-squares fit minimising S n i(P n i — P n ,l) 2 > 
for g modes of degrees I = 1, 2, 3, 4 and radial order n m ; n < n < n max , allows the values 
of Po> Vi, V2, and $, to be estimated It has been shown that by varying the values of 
n min an d fimax - for example n m ; n from 10 to 15 and n max from 27 to 35 - and by 
taking the mean values of the results for the global quantities, it is possible to deter- 
mine from a set of numerical periods the quantities Pq and V\ , related to the structure 
of the inner radiative zone. The values of Po and Vi are then shown to be very close 
to their asymptotic values. 

Similar good agreement is shown between the numerical periods P n j and the 
asymptotic periods P rh i- The fractional agreement is of order 10~ 3 or less for radial 
orders larger than 10, that is for frequencies less than 60 fiRz at I = 1, and for frequen- 
cies less than 100 /iHz at I = 2. At lower radial orders, the effects of the third-order 
terms and of the non-constant i9, which are neglected here, may account for fractional 
differences of up to 5 10 -3 . 

These asymptotic properties of g modes may be exploited in attempts to detect the 
modes, i.e., by searching for signatures of near-regular patterns in period, as we shall 
discuss in detail in Section 4. 



2.1.3 Effect of rotation : splittings, rotation kernels 

The rotation of the Sun induces a splitting of each multiplet frequency v n ; into 2Z+1 
separate singlet frequencies v n i m , where m is the azimuthal order of the oscillation 
(—1 < m < I). We define the splitting as: 

a n,l,m = { v n.l,m ~ v n,l,o)- 

The effects of the Coriolis force and distortions of the model must also be considered. 
The angular velocity fi is small compared with the characteristic frequency Og = 
^/GM/R 3 ((2/Q s ~ 5 x 10 3 ), so that the distortion of the model, which depends on 
(J?/J? g ) 2 , may usually be neglected. Generally a first-order analysis in the ratio Q/w is 
adequate to compute the rotational splittings (|Cowling fc Newing|ll949l : lLedouxlll95ll ). 
In this first-order approximation, the splittings in an inertial frame are related to the 
rotation by 

_ m J Q£- £* p&V - if (£x. £*) ■ fip&V 



(16) 



where the integrals are taken over the whole solar volume, dV being an element volume. 
Moreover, it is sufficient to adopt the corresponding (unperturbed) eigenfunctions £ 
of the nonrotating stellar model to evaluate the integrals. If the magnitude O of the 
angular velocity O is constant on spheres the splitting can be represented by the simple 
weighted averages: 

o~n,l,7n = m n (r)K.^' l (r)dr, (17) 
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Fig. 5 Rotational kernels /Cq' , caused by a rotation constant on spheres (cf. Eg. I17J I. plotted 
as a function of radius t/Rq. The quantity Rq x Kq' is plotted for modes 1 = 2 with radial 
order n = 5 (full line) and n = 1 (dashed line) in the upper panel and n = 20 (full line) and 
n = 10 (dashed line) in the lower panel. The dotte d lines represent the asympt otic envelope 
given by the function (P a /n 2 )(l - l/L 2 )N/r (from iBerthomieu fc Provostlfl993) . 



with 

Some typical examples of ICq' 1 are plotted in Figure [5] they have large amplitudes in 
the solar core, implying a high sensitivity to core rotation. 

As noted previously, low-frequency g modes of low degree have much larger dis- 
placements in the horizontal direction than in the radial direction, so that the sim- 
plified expression for the rotational kern els depends on l y on Using the asymptoti c 
expressions for the eigenfunctions (e.g. iTassoull Il98fj| ; Provost fc Berthomieul Il986l ). 
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Berthomieu fc Provostl ( 199lh have shown that the asymptotic kernels have a sensitiv- 
ity envelope in the radiative zone that is proportional to (1 — l/L 2 )N/r. This demon- 
strates the sensitivity of the frequency splittings of the gravity modes to the solar core 
rotation (Fig [5]). 

We know from the inversion of observed p-mode frequency splittings that the so- 
lar rotation is not dependent only on radius, but that the rotation in the convec- 
tive zone varies with l atitude as much as the rotation observed at the solar surface 
( Thompson et al.1119961 ) - 

If we take a simple law that reflects this differential rotation namely Q — Oo (r) + 
Qi (r)cos 2 (tf), we may then express the frequency splittings cr n i m of the g modes by 



O 1 



1 \ — 



JJ2 ) + fa I C W 



(19) 



with 



21 + 1 



{I + l)' 2 - m 2 
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where fi$ and Q\ are mean values of the rotation coefficients seen by g modes, and 

jQo,i(r){N/r)Ar 



(20) 



J?n_i = 



/ (N/r)dr 



(21) 



An important point to be made about the equations above is that <T n i m ~ mi?o(l — 
■j^), so that, at low frequencies, the splittings of I = 1 g modes are about half the size 
of the splittings of I = 2 g modes. This must be taken into account when attempts are 
made to identify g modes in the low- frequency range. When g modes with very low 
frequencies are considered, the ratio fl/uj is no longer small and another asymptotic 
analysis has to be made to take into account the fact that Coriolis effects cause the 
oscill ation eigenfunctions to be no longer described by a sing l e spherical harmonic 
(e.g. iBerthomieu et alJ 19781 : iDziembowski fc Kosovichevl 1 19871 : iDintrans fc Rieutordl 
|2000| ;1 Mathis et alj|2008l ). These modes can be expanded in terms of a hierarchy of 
unperturbed eigenfunctions £, a procedure which is valid because f2/O g is small even 
though fi/uj is not. However the usual first-order expression for the frequency splittings 
is rel evant for modes of frequency larger than about 20/iHz ( Berthomieu fc Provostl 

diil). 



2.2 g-mode frequency predictions and their reliability 

In order to aid observational attempts to detect and identify low-frequency p and 
g modes in the frequ ency spectrum, sev e ral groups have computed theoretical g- 
mode frequencies ( e.g. IProvost et all l2000l ; ICouvidat et al.ll2003l : llVIathur etalll2007l : 
ICox fc Guzikll2004l ). The reliability of these predicted frequencies depends first on the 
uncertainties of the solar modelling and on the sensitivity of the low-frequency oscil- 
lations to these uncertainties, and second on the numerical precision of the frequency 
calculations. Another important uncertainty lies in the estimation of the frequency 
splittings of the g modes, because these splittings depend of course on the rotation 
close to the centre of the Sun, which is not well constrained observationally. 
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2.2.1 Modern solar models and their uncertainties 



The large numb er of p-mode frequencies already detected lead to a seismic model 
of the Sun (e.g. iGoughl |2004 ) that is found to have a structure that is very close 
to the structures shown by modern solar models (e.g. IChristensen-Dalsgaard et al.l 
Il996l : iGoueh et alJll99d ). The only regions where the structure inferred by p-modc 
observations of the p modes remains uncertain is in the central few percent of the solar 
core. (Uncertainty in the angular velocity, as inferred from the frequency splittings of 
p modes extends through a large part of the solar core.) 

Standard solar models are assumed to be spherically symmetric, in hydrostatic 
equilibrium with no macroscopic motion except convection, no mass loss nor accretion, 
no rotation and no magnetic field. The structures are determined by the assumed 
microscopic properties of the solar material (equation of state, opacity, nuclear reaction 
rates etc.), and by the d escription of energy t ransport by convection usually (mixing- 
length prescription, e.g. Bohm- Vitensel [l958l : ICanuto fc Mazzitellilll99ll ; iGoueh et al.l 
19931). The num erical codes leading to these solar models have 



1996: lTurck-Chieze et al 



been carefully compared ( Monteiroll2009l ) . Solar models are calibrated by adjusting the 
initial abundances Y[ and Z\ of helium and heavy elements, and (usually only a single) 
mixing-length parameter which determines the efficacy of the convective transport in a 
mixing-length description, to obtain, for a one-solar-mass model, the solar radius and 
the solar luminosity at the solar age, and the obser ved surface metallicit y (Z/X)q, 
where X = 1 — Y — Z is the abundance of hydrogen. lAntia fc Basul (|2006h attempted 
to determine the solar heavy-element abundance without reference to spectroscopy by 
using helioseismic data. 

The introduction into the models of microscopic diffusion and settling of chemi- 
cal elements decreases the surface helium content, reducing the discrepancy between 
the observed and modelled sound-speed profiles, and also bringing the surface he- 
lium content, and the conv ection zone d e pth, of the models into better agreement 
with the observations, (e . g. lBurgerslll969l: ICox et al.lll989l; IProffitt fc Michaudlll99ll: 
Michaud fc Proflitt|[l99a ; IChristensen-Dalsgaard et al.lll993l : iThoul et al.lll994l ; llVIorell 



19971 ; iBrun et al.lll998l ). But what of the heavy-element abundances? 



The quantity of heavy elements assumed in the models is extremely important, 
because these elements dominate the radiative opacities. Usually one considers only 
the total heavy element abundance Zq or the so-called metallicity (Z/X)q. Until 
2004, the solar abunda nces were taken from resu lts given b y two rather close spectro - 
scopic analyses, namely [Grevesse fc Noelsl ( 19931 ) (GN) and lGrevesse fc Sauva 1 (1998) 
(GS). However, since then new heavy element abundances have been inferred using 
a 3D, NLTE analysis of the spectrosc opic data, which takes hyd rodynamics into ac- 
count in the transfer calculations as in lAsplund et ail (|2005l . |2009| . hereafter AGS) and 
ICaffau et all (|2009l ). 

Amongst other differences, the AGS abundances of C, N, O, Ne are substantially 
lower than the older GN and GS abundances, by 0.16 dex, 0.19 dex, 0.21 dex and 
0.24 dex. The reduction lowers the total surface metallicity of (Z/X)q (0.0166, com- 
pared to 0.0245 for GN and 0.023 for GS abundances). These changes have led to an 
extensive discussion in the literature on the robustness of the new results, which is 
beyond the scope of this paper, (e.g. Basu fc Antia| |2008l ; ICaffau et al.l 120081 ). In 



par- 



ticular, it led to an attempt by lAntia fc Basul ((2006J) to determine the heavy-element 
abundance seismologically without r eference to spectrosco py, by calibrating the average 
effect of ionization on 71. However. lAntia fc Basul ( 20061 ) simply adopted a particular 
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equa tion of state which, like all others, was known not to be exact (cf. iBaturin et al.l 
bood ). so their result cannot be relied upon. 

Solar models constructed with the GN or GS abundances in the generally accepted 
opacity formulae are in good agreement with helioseismic constraints. For example, 
the agreement for the sound-speed profile is within 0.2 perc ent, as shown in Figure 6 
(e.g. IChristensen-Dalsgaard et al.l [l996l ; iGough et al.|[l996h . In attempts to account 
for the discrepancies that remain, several non-standard physical processes have been 
included in the solar modelling, in particular processes to improve the description of 
the layers located just beneath the convection zone where the most marked differences 
between the observations and the models occur. We know from helioseismology that 
the transition from differential rotation in the convection zone to almost uniform ro- 
tation in the radiativ e zone below occurs in a shallow layer called the tachocline (e.g . 



ISpiegel fc Zahnlll992l ; iKosovichevI Il996l ; iThompson et al.lll996l ; | Corbard et al 



The shear of the rotational velocity necessarily induces mixing (I Hughes et al 



1999). 



2005). 



Various descriptions of this mixing have been included in some models, leading to re 
ductions in the dif f erences between the observations and models, (e.g 



Richard et alj 



199rj lGabriel|[l997l ; I Morel et al.lll99Sl ; iBrun et al.lfl999l ; lElliott fc Gougbll 19991 ) 

Other non-standard processes, such as overshoot of convective elements at the 
boundary of the convective core which appears at the end of the pre-main-sequence 
evolution, or mass-loss occurring during the initial stages of solar evolution, have been 
studied, with special emphasis placed on the changes induced in the most central parts 
of the models which are not well constrained by observations of p modes (see below). 

The AGS abundances have caused considerable concern. This issue is simple to 
state: Except in the near-surface layers, which are of no import to this discussion, we 
already know from seismology the sound speed and the density stratification through- 
out the Sun to a precision easily good enough to be confident in the inference. If 
one then accepts the basic tenets of solar evolution theory - namely that a helium- 
abundance augmentation has occurred in the core due to nuclear transmutations on 
the main sequence and that the Sun is in thermal balance, energy being transported 
throughout the entire region beneath the convection zone by radiative transfer - and if 
one accepts also the equation of state (which we surely know to well within the limits 
set by the abundance change) and the nuclear reaction, then one can infer the absolute 
value of the hydrogen-helium abundance ratio and hence the temperature stratification 
(by demanding that the energy-generation rate in the core is equal to the luminosity 
observed at the surf ace) from wh ich can be inferred the value of the opacity beneath 
the convection zone |Goughll2004l ). The issue, therefore, is simply how to reconcile that 
value with the photospheric abundances. At first sight, the apparently most straightfor- 
ward resolution is that the abundances in the convection zone differ from these in the 
radiative interior, presumably as a results of accretion of metal-deficient material after 
the Sun has (almost) reached the main sequence and had a relatively shallow convec- 
tion zone, which is co ntrary to standard a ssumptions. We should point out, however, 
that iGuzikl (|2006l ) and lCastro et al.1 (|2007l ) have failed to construct such a model that 
is in adequate agreement with helioseismology. 

It follows that using the AGS abundances in the accepted opacity formulae leads to 
solar models which disagr ee with the helioseismically determined constraints on the so- 



solar models wmcn disagree witti tnc lielioseismically determined constraints on trie so- 
lar interior properties (e.g. Turck-Chieze et al. 2004; Montalban et al .. 2004; Ba hcall et al.l 



2005; iGuzik et al.l 120051 ; IZaatri et al.l 120061 . 120071 ). There is a large discrepancy in the 
sound speed (e.g. Figure|6}. The surface helium content of these new models is too low, 
and the convection zone is too thin, compared to the observed (helioseismic) values. 
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Fig. 6 Relative density profile and relative squared sound-speed di fferences betwe e n stan dard 
models using a new estimate of the CNO composition proposed bv lAsplund et al ] p005h (full 
line with seismic error bar) and the heliosei smic observation s. The dot-dash line corresponds 
to a model computed with abundances from IrTolwcgcr (2001]), intermediate between AGS and 
GN. The seismic model (thin line) is an model which reproduces properly the observed sound 
speed and is used to predict gravity modes and neutrino fluxes (After Turck-Chiczc ct al. 
12004 ). 



Many other attempts, and suggestions, h ave been ma de to try to reconcile the AGS 
abundances with the helioseismic results (e.g. lGuzfklboOfjl). Evidently, increasing the ra- 
diative opacities below th e convection zone (e.g. ISerenelli et al ] |2004l ; IMontalban et all 
|2004| ; iBahcall et ail 12004" ) decreases the discrepancies if one insists on lowering the 
heavy-element abundances to t he AGS values. Using artificially enhanced diffusion 
rates for He and heavy elements ( Basu fc Antial2004l ; IMontalban et al.ll20o3 ; lGuzik et al] 
l2005l ; lYang fc Bill2007l ). does improve the agreement with seismic sound speed, though 
such changes by the amount required have not been justified. Rotational mixing, acting 
in the tachocline just beneath the convection zone can lead to models with acceptable 
valu es of the surface helium abundance, but their convection zones are too shallow 
(e.g. iTurck-Chieze et aljliooj ). Taking account of tachocline mixing (perhaps calling 
it convective overshoot) helps to bring the convection zone depth closer into line, but 
the discrepancies shown in the sound-speed profile are not necessarily reduced in size 
(e.g. IMontalban et al.ll2006l ). 

In sum, at the time of writin g the solar abundan ces remain a major source of 
uncertainty in solar modelling (e.g jBasu fc AntialboO&l and references therein). 

Finally in this section we consider mechanisms to describe transport of angular mo- 
mentum and chemical elements which are not fully accounted for in solar and stellar 
evolution codes. The respective roles in angula r momentum transport of dynamical pro - 
cesses involving rotation, mag netic fields (e.g. iGough fc Mcintvrelll998l ; lGaraudll2002r) 
and gravito-inertial waves have been much discussed in the solar conte xt ( Schatzmanl 
Il993l . Il99fj ; iTalon fc Zahnlll99Sl ; iTalon et al.l |2002| ; iMathis et al1l2008l ; iGoughl l2009f T 
When only the meridional circulation and the "classical" hydrodynamic instabilities 
are invoked, models predict a Sun with a radiative interior and core rotating much 
faster than the surface (e. g. IPin sonneault et al.l ll989l). which is incomp atible with p- 
mode helioseismology (e.g. lThompson et al.lll996l ). lMathis fc Zahnl (|2005l ) have recently 
derived, in a self-consistent manner, a formalism to describe the effects of an axisym- 
metric magnetic field on the meridional circulation and the turbulence generated by 
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Fig. 7 Differences 8u = i> n ; — u n ; between eigenfrequencies f n ; computed by an oscillations 
code, and predictions of the frequencies given by the integral expressions for a reference solar 
model, p n j [cf. Eq. I|10j l] for modes of degree Z=0 (dots), 1 (triangles), 2 (full stars) and 3 
(squares). 



the shear of differential rotation. This formalism is now implemented in some evolution 
codes. 



2.2.2 Numerical precision 

The numerical precision of the frequencies depends both on the physical and numerical 
accuracy of the evolutionary models and on the algorithms used to compute the oscil- 
lation frequencies themselves. Many evolution codes have been compared, in the 1990s 
for the Sun (GONG solar model comparison) and more recently in a stellar context 
as part of the CoRoT Evolution and Seismic Tools Activity (ESTA) ( Lebreton et al.l 



l200al : lMonteiroll2009l ). The numerical precision of the frequencies may be estimated 
for a given stellar model by considering the range of v alues obtained fro m different 
oscillations codes, which compute the mode frequencies ( Mova et al.ll2008h . 



For a given oscillations code, an internal accuracy may also be estimated by com- 
pari ng a numerical eigenfrequency v n ; with its integral value, as obtained from Eq. (|10|l 
(e.g. IProvostJliooil ). Any inconsistency, either in the computation of the oscillation fre- 
quencies or in the equilibrium models, gives rise to a non-zero value of v n i — v n i as, 
of course, does the neglect of the surface term in Eq. (|10[) which are actually negligible 
only at low frequency. 

For g-mode calculations one also needs a solar model with a sufficiently large num- 
ber of mesh points (where the model quantities are given explicitly), well distributed 
in solar radius, i n particular close to the centre of the model where g modes have large 
amplitudes (e.g. IChristensen-Dalsgaard fc Berthomieulll99ll ). Figure [7] shows that, for 
a reference solar model, a numerical precision of g-mode frequencies - as measured 
by v n i — v n i - better than 0.01 /iHz for frequencies lower than 200 /iHz has been 
achieved. This difference is much larger than uncertainties we expect to obtain from 
observations (should the g modes be detected). Therefore the numerical procedures 
should be improved. However, as we shall show below, this error is much smaller than 
errors introduced by uncertainties in the theoretical structure and the dynamics of the 
models. Therefore g modes have diagnostic potential. 
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Table 1 Asymptotic coefficients of g modes, i.e., Po, Vi, V2, and Plotted are mean values 
of the coefficien ts for n m j n from 10 to 15 and n max from 27 to 35, as computed f or models 
Ml to M8 from lProvost et al.l (|2000|> (1), M-GN & M-AGS from lZaatri et al.l J2007I) (2). The 
estimation of the uncertainty, in percent, was obtained from the scatter of the least-squares fit 
values around their mean values. In this table the following differences show the following effects 
(for details sec Provost et al. (2000)): (M6-M1) shows the effect of increasing age from 4.65 to 
4.75 Gyr; (M2-M1) shows the effect of increasing the surface metallicity from 0.0245 to 0.0260; 
(M3-M1) and (M8-M1) show the effects of mass loss, set to a rate of M = -5 X 10~ 10 Mq?/~ 1 , 
occurring for 200 Myr at the beginning of the e volution, with additional turbulen t mixing below 
the convection zone in model M8, according to Montalban & Schatzman (2000); (M4-M1) and 
(M5-M1) show the effects of overshoot of, respectively, 0.2 and 0.25 pressure scale heights at 
the edge of the convective core which appears at the end of the pre-main sequence evolution; 
(M1-M7) shows the effects of improved nuclear reaction rates; (M-AGS - M-GN) shows the 
effects of changing the chemical abundances from GN to AGS (see text); and finally (M-AGS 
- M-GN) shows the effects of changing the chemical abundances from GN to AGS (see text). 





Po 

(mn) 


% 


Vi 


% 


v 2 


% 




% 




Ml 


35.08 


0.04 


0.437 


1.6 


5.69 


4 


-0.160 


11 


(1) 


M2 


34.94 


0.04 


0.439 


1.6 


5.67 


1 


-0.161 


11 




M3 


34.63 


0.04 


0.444 


1.6 


5.72 


1 


-0.161 


11 




All 


35.39 


0.04 


0.430 


1.6 


5.64 


1 


-0.161 


11 




M5 


36.75 


0.03 


0.399 


1.1 


5.46 


3 


-0.164 


8 




M6 


34.77 


0.04 


0.440 


1.6 


5.71 


1 


-0.161 


11 




M7 


35.42 


0.04 


0.433 


1.6 


5.61 


1 


-0.162 


11 




M8 


34.78 


0.07 


0.437 


1.3 


6.12 


6 


-0.217 


14 




M-GN 


35.09 


0.03 


0.442 


1.5 


5.77 


3 


-0.166 


9 


(2) 


M-AGS 


35.69 


0.04 


0.433 


1.4 


5.96 


3 


-0.157 


11 





2.2.3 Sensitivity of low-frequency g- and mixed mode frequencies to solar models 

As a first step, the sensitivity of g-mode frequencies to uncertainties in the solar models 
may be studied by taking advantage of the asymptotic properties of g modes. As noted 
previously, the three global coefficients Po, V\, V2, and the phase factor i9, which 
characterize the frequency may be obtained by a least-squares fit to a set of numerical 
frequencies. The values of Po and V\ are typically very close to those derived directly 
from the integrals (cf. Eq. (|11[) ). and results for some models are given in Table [1] 
In this table one can see the effects of changing some fundamental model parameters, 
such as the age to which the model is evolved, the metallicity, and also some important 
physics such as the nuclear reaction rates. 

We find that the values of the asymptotic coefficients are very similar for models 
with slightly different fundamental parameters and different physics, including changes 
between the ol d and new metall i cities. Ignoring gravit ational settling increases Pq by 



between the old and new metallicities. Ignoring gravitational settling increases Po by 
about 0.7 min jMorel et al.lll997t ). lMatmrr et al]||2007l ) have found that the theoretical 



value of Po changes by no more than 1 min when models with plausible differences 
in treatments of microscopic diffusion, turbulence in the tachocline, and the chemical 
composition, are considered. The results give an idea of the possible range of values for 
these asymptotic coefficients. 

The sensitivity of the frequencies to uncertainties in the models can also be studied 
in more detail by looking at t he relative frequency differences compute d for two differ- 
ent evolutionary models (e.g. IProvost et aUbOQOl ; iMathur et a.1.1 l2007h . Some general 
properties can be seen in Figures [8] [9] and 1101 (i) at low frequencies, the frequency dif- 
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Fig. 8 Relative frequency differences, 5u/u, between a solar model and the reference model 
Ml. Results are shown for g modes (open symbols) and low-frequency f and p modes (full 
symbols) with degrees I = 1 to 6, to test the sensitivity to the sola r parameters: age (m odel 
M6 - St = O.lGy), and luminosity (S L/L=3.6 10 -3 ). Adapted from lProvost et al.l (120001') ■ 



ferences between two models obey approximately the asymptotic relation of Eq. (|14l) ; 
(ii) the largest frequency differences appear in the range 200 to 300 /iHz, where there 
are mixed modes, and the behaviour is more irregular as a function of frequency than 
the behaviour shown at very low frequencies; and (iii) in contrast, the sensitivity of the 
low-frequency p-mode frequencies to changes in the solar models is much smaller than 
for the g modes. 

The sensitivity of low-degree g-mode frequencies to changes in the age and lumi- 
nosity of solar models is shown in Figure [5] The frequencies are seen to increase with 
increasing age, increasing metallicity, Z/X, and increasing luminosity. For instance, an 
increase of either 1 % in age, 6 % in metallicity increases the frequencies by about 0.4%. 
An increase in lumino sity of 3.5 %, which corresponds to the upper value reported in 
iGuenther et al.l (|l992h . increases the g-mode frequencies by 0.2%. A frequency augmen- 
tation of abou t 1 % arises from augmenti ng the nuclear reaction rate s above the values 
in the tables of lCaughlan fc Fowled ( 19881 ) to the most recent tables of lAdelberger et alj 
( 19981 ). 

A reduction in the solar abundances from the old (GN/GS) to the new (AGS) 
values reduces the frequencies of the or der of 1.5% or more, as shown in see Figure [5] 
( Mathur et~aH 120071 ; IZaatri et al.l l2007h . This effect is in agreement with the global 
asymptotic p arameter resu l ts sh own in Table [T] for models M-AGS and M-GN. We 
also note that IZaatri et all ( 2007) have reported that the characteristic low- frequency 
g-mode period Pq decreases by increasing the neon abundance i n the new solar mixtu re 
(Figure |§J) which is consistent with inferences from the work of lMathur et al.l ( 20071 ). 

The effects of differe nt descriptions of mixing in the tachocline have been considered 
bv lMathur et all d2007l). who showed that the introduction of horizontal diffusion, as 
suggested bv lSpiegel fc Zahnl (| 19921 ). would decrease the g-mode frequencies by about 
0.25%. Core overshooting and mass loss during initial stages of evolution decreases 
and increases the g-mode frequencies, respectively, by about 1 % (Figure fTU|) . 

In conclusion, the sensitivity of the g-mode frequencies to uncertainties in the solar 
models appears to be of order at least 1% for the cases considered here, which corre- 
sponds to frequency differences of the order 1 /iHz around 100 ^iHz and 2 /iHz around 
200 ^tHz, i.e. much larger than the numerical precision of the frequency computations 
(down to 100 ^iHz). 
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Fig. 9 Relative frequency differences between the g modes and the gravest low- frequency 
p modes of a calibrated solar model computed with the old GN abundances (M-GN) and a 
model computed with the AGS abundances (o). The other symbols, (A, *, □, X, •) correspond 
to models computed with increasing changes to the the neon abundance of the models with a 
dex of 8.10, 8.29, 8.35, 8.47 and 8.35, respectively. From lZaatri et al.l J2007f) . 
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Fig. 10 Relative frequency differences, Sv/v, between a solar model and the reference model 
Ml, shown to test the sensitivity to the following physics (see caption of Table 1): overshoot of 
0.2 pressure scale heights (mo del M4); and s t rong mass loss rates (model M3) (same symbols 
as in Figure [§). Adapted from lProvost et all (120001 ') 
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Fig. 11 Splittings of low-degree, low-frequency sectoral m = I modes as a function of the 
frequency for a reference model and for a simplified rotation law, which mimics for r > 0.2 Rq 
the solar rotation inferred by inversion of the observed p-mode frequency splittings. Shown are 
predictions for different values of the rotation i? c below 0.2 Rq. Left panel: Q c = 0.433 ^Hz. 
Right panel: Q = 2Q C for r < 0.2Rq (open symbols: g modes; full symbols: f and p modes). 
From lProvost et al.l d2000T l. 



2.2.4 Sensitivity to the rotation and dynamics in the core 



The sensitivity of g-mode frequency splitting to the dynamics of the solar core may be 
investigated by the forward method, using a solar rotation law, which, for r > 0.2 Rq, 
mimics the rotation rate obtained by inversion of the p-mode frequency splittings, 
i.e. rigid rotation below the convection z one and latitudinal differential rotation inside 
the zone (see e.g. ICorbard et ai1ll997l ; iDikpati et al.ll2002h . Since the core rotation 
is very uncertain, analyse s of this type assume a constant rotatio n in the core r < 
0.2 710 equal to O c (e.g. IProvost et al.l [2OO0I ; iMathur et al.l 120081 1 . Figure ITT1 shows 
the variation with the frequency of the model-predicted splittings, o~ n n, of sectoral 
g modes, for two different values of Q c : either the core rotation is assumed the same 
as in the radiative zone (J? c = 0.433pHz; left panel), or it is assumed to be twice this 
value (right panel). At low frequencies, the splittings depend on the degree I, with a 
behaviour close to the asymptotic behaviour of g modes, i.e. a splitting proportional 
to 1 — 1 / L 2 . At high frequencies the splittings instead tend to the asymptotic p-mode 
values, with a very weak dependence on the degree I arising from differential rotation 
in the convection zone. As expected, if the core rotation is increased by a factor of 
two, the g-mode splittings are increased significantly by a factor of about 1.75, while 
the p-mode splittings are almost unchanged. The splittings of the mixed modes have a 
complicated dependence on degree, frequency and the shape of the rotation law below 
0.2 Rq. If the rotation of the solar core differs notably from that of the radiative zone, 
such behaviour will make their detection and identification in the observed spectrum 
harder. 

The internal magnetic field, which is often neglected, can also split the frequen- 
cies into se veral components and thus it can affect the reliability of predicted g-mode 
frequencies. iRashba et al. I |2007l ) attempt to estimate analytically the effect of a mag- 



netic field in th e radiative zone o n the g-mode frequencies, using a MHD perturbative 
correction (e.g. lUnno et al.lll989l ). They show that a 1-per-cent g-mode frequency shift 
can be obtained for magnetic field as low as 300kG in the radiative zone, for modes 
of radial order n=20 and degree 1=1 (y ~ 33/iHz), and they argue that a similar shift 
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Fig. 12 Equatorial rotation profiles below 0.4 Rq reconstructed by inversion of artificial data, 
using eight g modes, each with an assumed frequency splitting uncertainty of 7.5 nHz. Results 
are shown for three different rotation laws in the solar core (dashed line rigid, full line step 
profile, long dashed line smooth profile). The continuous blue, red and dotted li nes represent, 
respec tively, the input step, smooth and rigid artificial rotation profiles. From lMathur et alj 
J2008D . 



for the case of low-order g-mode frequency cannot result from a central magnetic field 
unless its strength exceeds 8 MG which would be a very large magnetic field. 

The rotational kernels would help to infer the core rotation by inversion, once the 
splitting of g modes could be measured. Preliminary work using inverse methods has 
been carried out to try to know how measured splitti ngs of some g modes would influ- 
ence o ur knowledge of the rotation of the solar core (jMathur et al.ll2007l ; I Garcia et al.l 
l2008rJ ). Numerical inversions using artificial data, corresponding to different laws of 
rotation in the solar core (see Figure I12[) were performed. It was shown that the in- 
troduction of only one g mode gives a much better idea of the rotation profile. The 
introduction of eight g-mode frequencies with an error bar of 7.5 nHz in the data allows 
one to recover the rotation profile in the core quite well, with the low-degree p modes 
contributing to constrain the solar rotation profile down to 0.2 Rq. As expected, with- 
out g modes, it is not possible to distinguish between the three assumed core rotation 
laws 



2.3 Conclusion 

It appears that g-mode frequency predictions are affected by two sorts of uncertainty. 
First, the faults in the solar modelling induce uncertainties in the solar stratification, 
hence in the Brunt- Vaisala frequency which determines mainly the g-mode frequen- 
cies. Some contributions have been estimated to 1% (see subsection 2.2.3), but the 
consequences of the lack of a reliable description of the dynamical and magnetic effects 
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remains to be estimated. Second, the prediction of the splittings depends on the core 
rotation, which remains poorly known. A rotation twice larger below O.2i?0 would 
increase the low-frequency g-mode splittings by a factor about 1.75. 



3 Excitation and damping of g modes 



In the past, s everal attempts have been made to estimate solar g-mode ampl itudes 
( Gough|[l985l : iBerthomieu fc Provost] Il990l ; lAndersenl fl99tjl : iKumar et al.lll99rJ ). The 
common res ult of these prediction s is that g-mode amplitudes are not larger than a 
few mms -1 jElsworth et aLlbood ). 

The investigations of g-mode amplitudes first focused on the linear stability of the 
modes so as to determine whether or not they could be excited by overstability (i.e. 
thermal instabilities). Results published to date suggest that g modes are most likely 
linearly stable (Section l3.2p . Consequently most of the theoretical efforts have therefore 
assumed that g modes are intrinsically damped and excited stochastically by turbulent 
convection (Section 13. 4[) . as it is also believed to be the case for p m odes, whether 
the excitation is predo mina ntly in the body of the convection zone jGoughl Irjjj 
IBerthomieu fc Provostlll990l; IKumar et al.lll996l: iBelkacem et al.ll2009l ). or in the lower 
overshooting layers (e.g., lAndersenl 1 19961 : iDintrans et a.1.1 120050 . Note also that other 
mechanisms such as mode coupling (e.g., Dziembowski 1983 HGuenther fc Demarquel 



mechanisms such as mode coupling (e.g., Dziembowski laso; (jucnthcr &i Dcmarquc 
1 19841 : lAnddll98d:IWentzeJll98^:IWolff fc O'Donovanll2007r ). or excitation by magnetic 
torques (|Dziembowski et al.lll985T ) have been investigated. 

The resulting amplitudes are determined by a balance between the intrinsic linear 
stability of the modes, which are characterised by the damping rates j), and the energy 
input via nonlinear interactions with the turbulence. These will be discussed separately 
in Section I5~2l and I5~4l 



3.1 The separation of the oscillation modes from the convection 

Before describing some of the details of what has appeared in the astrophysical litera- 
ture, it is instructive first to ask just what the question is that we wish to ask. The glib 
response is simply: to what amplitudes does the turbulence drive the modes? But before 
we can address that matter, we should try to be clear in our minds how to separate the 
motion into convection and modes of oscillation. There have been several attempts to 
accomplish that task, mainly unpublished because unfortunately unambiguous success 
has never been achieved. Therefore we shall be brief. 

In the case of radial (p) modes, superficially the issue is reasonably clear: the modes 
are horizontally uniform, and convection, which is motion driven by buoyancy, is intrin- 
sically horizontally nonuniform. Therefore one simply takes a horizontal average, and 
regards the average motion as the oscillation mode and the remainder as convection. 
This is a fairly s afe procedur e, at least for isolating radial modes (and was justified to 
some degree by lGougldll969h . provided that one takes the horizontal average of pw, 
where p is density and w is the vertical (radial) Eulerian component of the convec- 
tive velocity, rather than w, to define the oscillation mode. Incidentall y, convection , 
even w h en approxim ated by a mixing-length-type theory of eddies (e.g. lGough|[l965l . 
Il977al lbl; lunnolll967n , is most naturally described in (local) Eulerian coordinates; since 
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the eddies maintain a degree of integrity, they are advected by the large-scale pulsa- 
tional flow, which is therefore most appropriately described in Lagrangian coordinates. 

Separating nonradial modes from convection is a more difficult task. It can be 
accomplished only approximately, and that only when there is a good separation of 
scales. Modes of low degree have horizontal scales comparable with the radial coor- 
dinate r, and because most of the energy in the convective motion is in scales rather 
smaller than that (at least on timescales comparable with the periods of the p modes 
and grave (n < 3) g modes) one can at least separate the motion, this time by taking 
appropriate large-scale spatial and temporal a verages. In this way, a mixed Eulerian - 
Lagrangian coordinate system can be set up jGabriel et al.1 1 19751; llJnno et al.|[l989l ) 

"a 



similar to that defined originally for radial pulsations 1 Goughl 19691 ). However, such a 



separation cannot be achieved for the very-high-order g modes that resonate with the 
turbulence at the base of the convection zone. In that case it is n ecessar y to re sort 
to more involved procedures. For example, one may try, following IPovetl (| 19831 ). to 
decompose the motion into a superposition of linear p modes and g modes, regard- 
ing the function space spanned by the direct, so-called g + , modes as convection; that 
procedure suffers from ambiguity when trying to estimate the overlap integrals that 
couple the two kinds of motion, because there can remain some freedom in how the 
decomposition is carried out. Most of the analytically based studies simply separate 
the motion by fiat. And now, having issued our warning, we shall do likewise. 

Before proceeding, a word concerning the analysis of numerical simulations is in or- 
der. Simulations have been carried out to shed light directly on the mode-excitation is- 
sue, although in some cases al so with the intention of calibrating analytical approxima- 

tions to the energy-input rate (jRosentha lll998l ; IStein fc Nordlundl200ll ; lRogers fc Glatzmaierl 
l2005l ; iRogers et alJuOOa ). Because convective motion cannot be unambiguously sepa- 
rated from what we imagine to be the oscillation modes, the details of what one envis- 
ages to be energy conversion cannot be isolated. But what one can do is to look in the 
far field where convection is ignorable, provided the domain of the simulation is exten- 
sive enough, and there determine the spectrum of the g-wave and p- wave radiation by 
projection onto linearized wave functions (e.g.. iDintrans et al.ll2005l ). 



3.2 Linear g-mode stability 

The d iscussion of the stability of g modes in a star has a long history. [Ledoux fc Sauvenier-Goffin] 
( 195d ) were probably the first to address the general problem, although they limited 
their discussion to white dwarfs. As is the case for acoustic modes, the stability of 
gravity modes depends on the integrated effect of various physical mechanisms which 
can either drive or damp the oscillations. The excitation and damping processes usually 
take place in layers with rather small radial extent: these are the ionisation regions of 
hydrogen and helium where the Eddington valve can operate, the energy-generating 
core where temperature-sensitive nuclear reactions can feed energy into the oscillatory 
motion, and also the highly superadiabatic upper convective boundary layer and, for g 
modes, the vicinity of the boundaries of convection zones where dynamical interaction 
between the convection and the oscillations is relatively strong. In the following section 
we review these mechanisms. 
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3.2.1 Linearized equations governing g-mode oscillations 

One first takes a statistical (temporal) average of the governing equations (which deter- 
mines the basic, background state of the star) and subtracts it from the full equations 
to yield the fluctuation equations. Having separated the fluctuating variables into an 
oscillatory mode contribution and convection, one linearizes the fluctuation equations 
in the mode variables, and then typically arranges the outcome with the mode variables 
on the left and the fluctuating terms associated with the convection on the right. In 
all studies to date the mode amplitudes are regarded as being so small that they do 
not influence the inhomogeneous convective driving term on the right. Therefore each 
mode can be treated separately. The equation of motion can then be written formally 
as 

Cd = ${v,T'), (22) 

where £ is the displacement eigenfunction of the mode, with (v,T') is the convective 
velocity and temperature fluctuation, and complex frequency 2-KV n i m = uj — in (See 
Eq.[T}, in which u) is real and positive, n is real. The term 5 represents inhomogeneous 
stochastic driving and damping terms that depend only on convective quantities (i.e. v 
and T 1 ) unpertub ed by the oscillations. The spatial differential wave operator C (e.g., 
lUnno et al.lll989l for details) depends upon ui — \r\ and the background state of the 
star, which itself depends on the convective fluxes of heat and momentum; it depends 
also on the perturbation to those fluxes that is produced by the oscillations. Ignoring 
the right-hand side of Eq. (|22[) yields the equation of free oscillation (i.e. free from the 
inhomogeneous damping and driving, but including in principle the linearized momen- 
tum and convective momentum flux perturbations that are induced by the oscillations) , 
which in the adiabatic Cowling approximation can be approximated by Eq. (j4]). 



3.2.2 Excitation by the e mechanism 

The first dynamical investigation of g-mode instability of the Sun in which the nuclear 
reactio ns were perturbe d consi stently was carried out in the quasiadiabatic approxima- 
tion by iDilke fc Gougbl |l972f ) who approximated the eigenfun ctions in the radiative 



interio r by those of a polytrope of index 3. Soon afterwards iRosenbluth fc Bahcalll 
(1973) discussed the thermal instability of hydrostatic disturbances. The aim in both 



papers was to question the assumptions upon which standard solar models are built 
with a view to addressing the still-existing solar neutrino problem. Dilke & Gough 
found that low-order g modes could have been dynamically unstable in the early evo- 
lutionary stages of the Sun owing principally to the strong temperature dependence 
of particularly the 3 He( 3 He,2p) 4 He reaction which manifests itself in the energy bud- 
get when the reactions of the p-p chain are thrown out of balance by the dynamical 
oscillations, and also to the fact that the Sun has a convective envelope deep enough 
to provide an effective evanescent shield; earlier investigations confined to stars with 
radiative envelopes (e.g. lAure1ll97ll ) had indicated that g-mode amplitudes are so high 
in the surface layers that radiative damping would overwhelm any driving in the core. 
Rosenbluth & Bahcall found that their solar models are stable to non-radial thermal 
instability. 
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Fig. 13 Theoretical growth rates 77 of three low-order g modes as a function of age t in the 
early stages of the main-sequence evolution of the Sun. The u nits are 10~ 7 y —1 for r\ and 
the a ge is measured in units of 10 s y of the equilibrium models ( Christcnscn-Dalsgaard et al, 
1974). The current age of the Sun is about 4 10 years. Modes are found linearly unstable for 
77 < and stable for 77 > 0. 



3.2.3 Radiative damping 



Foll owinglDilke fc Goughl (Il972l ). lchristensen-Darsgaard et alJ l|l974h . lBourv et al.1 (|l975l ) , 



and iShibahashi et al 



uEh ( iy 

(|1975T ) addressed the g-mode instability using eigenfunctions of 
solar models, although still in the quasi-adiabatic approximation. All authors found 
unstable low-order g modes in models younger than the present Sun. The instability 
arose because the temperature sensitivity of nuclear reactions perturbed on a g-mode 
timescale is substantially greater than it is on a ste llar-evolution timescales - the latte r 
is too weak to overcome radiative damping (e.g. iDziembowski &i Sienkiewicd Il973l ). 
For the present Sun, however, all authors found all g modes to be stable, at least if 
the evolutionary consequence of earlier instability could be ignored. It is particularly 
so in the early stages of evolution when the star was les s centrally conde nsed that 
the g-mode amplitudes are relatively large in the core (e.g. iNoels et al.lll974h , and are 
therefore more receptive to nuclear driving. Figure [13] shows low-order g-mode growth 
rates for an evolving 1 Mq star computed in the quasi-adiabatic approximation. 

The first f ully nonadiabatic stability compu t ations of solar p and g modes were 
undertaken by IChristensen-Dalsgaard fc Goughl ( 1975h . The authors concluded that 



nonadiabatic processes in the very outer layers of the present Sun dominate over the 
destabilising influence of the temperature-dependent 3 He( 3 He,2p) 4 He reactions in the 
core, suggesting that low-order g modes are very likely stable. An important deficiency 
in these computations, however, is that the modulation of the convective fluxes of 
heat and momentum were omitted. iKosovichev fc Severnvl ( 19851 ) investigated the g- 
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mode instability for solar models with different metallicities and concluded that the 
instability is stronger in low metallicity models. 



3.2.4 The k mechanism 



In stellar layers where partial ionisation of chemical elements takes place, opacity can 
increase with temperature during the compression phase of a stellar pulsation, tending 
to produce an Eddington valve that can lead to instability. The radiative (energy) 
flux through such a layer is then absorbed by the stellar matter more effectively than 
it would have been otherwise. After the subsequent expansion phase the previously 
absorbed excess energy is released, the net effect over the whole cycle being to feed 
energy into the pulsation. This mechanism is respon sible for pulsations of some variable 
stars such as 8 Scuti, 3 Ceph ei and Cepheids (e.g. iDziembowski fc Pamiatnvkhlll993l : 
iGautschv fc Saiolll995l . fl99g ). This process is effective only if much of the heat flux in 
the relevant layer is carried by radiation, which in the Sun is not the case. However, 
the valve can operate also as a result of a perturbation to the convective heat flux. 
Irrespective of whether the flux is predominantly radiative or convective, in any star 
the valve is most effective if the thermal relaxation time of the material in and above 
the relevant layer is similar to that of the period of the pulsation. Because the thermal 
relaxation time increases with depth, low-period oscillation modes are driven in deeper 
layer s than are corresponding short-period oscillations (e.g. IPamvatnvkhll 19991 : IHoudekl 
l2003h . 



3. 2. 5 The effect of convection 



bai 3 (|l98Ch reported the first non-adiabatic solar g-mode stability computations that 



incorporated the mo dulation of the convection, adopting the treatment of convection by 
iGabriel et all [|l975h . which is based on Unno's (1967) time-dependent mixing-length 



formalism. In these computations, only the heat flux was taken into account: the mo- 
mentum fluxes (Reynolds stresses) were omitted from both the equilibrium structure 
an d the s tabili ty analysis. 

I Said ( 198(1 ) confirmed the earlier findings by IChristensen-Dalsgaard et al.l ( 19741 ) 



of g-mode instability during the early evolutionary phases of the Sun. Interestingly, 
for solar models with an age of 4.5 x 10 y he found that some of the g modes are 
additionally driven by the hyd rogen ionization zone in the very outer layers of the 
star. However, as discussed by iDilke fc Goughl |l972j), the main agent driving grave 
low-degree g modes appears to be the 3 He - 3 He reaction in the energy-generating 



It was reported bv lGoughl (|l965l . Il980h that the dynamical effects arising from the 



turbulent momentum flux perturbations contribute significantly to the growth rates 
of solar oscillations. Stability computations of a star in which both the convective 
heat flux and the Reynolds stresses were included, in both the equilibrium structure 
and th e pulsation ca lculations, using t he time-dependent mix ing-length formalism by 



lGoughl(|l965l . ll977al '). were reported bv lBaker fc Goughl (|l979T l for acoustic osc i llations 



in RR Lyrae stars. A subsequent detailed analysis of the Sun by iBalmforthl (1992b) 



revealed the manner in which the phases of the fluctuating Reynolds stresses and the 
density perturbations differ, causing r adial p modes to be damped. Similar findings 
were reported also by iHoudek et alj ( 19991 ) in other main-sequence stars. Different 



results were obtained from the computation of iDupret et al. (2006a). which suggest 
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that the perturbation of the convective heat flux i s dominant. The formalism that 
wa s adopted is based on th e approach of lUnnol Jl967) g e neraliz ed to non-radial modes 



bv lGrieahcene et al l rtiooiih . In addition, iDupret et al.l (2006a) found that the effect 



of the turbulent pressure modulation is partly compensated by the perturbation of 
the d issipation rate of turbulent kinetic energy into heat (see also iLedoux fc Walravenl 

GUI). 

For high-order g modes [y < lOOitHzl. iBelkacem et al" ll2009h have sh own, using 



a partially non-adiabatic non-radial computation ( Grigahcene et all |2005h , that the 
damping rates are dominated by radiative losses and are insensitive to t he convection- 
pulsat ion interaction. This result is in agreement with previous findings of lKumar et al] 
(|l996h . 

It is perhaps useful to point out the fundamental approximations that have been 
used to model the temporal modulation of the convection by the oscillations. Formu- 
lae for these modulations are needed when estimating the influence of convection on 
the frequencies and, more importantly, the growth rates of the oscillations. We limit 
ourselves here to those procedures that have been used for estimating solar g-mode 
amplitudes. 

The most naive procedure is to presume that the convective fluxes simply relax 
exponentially on a timescale r c towards the time-independent formula evaluated for 
the current (time-varying) environment. Thus, if F c is a component of a flux, be it of 
heat or momentum, it is taken to obey an equation of the type 

dF c _ F c0 - F c 

where r c is a multiple of i/w, I being the mixing length and w a characteristic convec- 
tive velocity, both of which are evaluated typically for the non-oscillating background 
state. The function F c q is some chosen formula for the corresponding flux evaluated 
for the instantaneous state of the star as though that state were not pulsating. Ev- 
idently, there is some formal freedom in the choice of that formula, because it was 
derived under the assumption that the background state is in (statistical) hydrostatic 
equilibrium, whereas when the star is oscillating that is not the case. Therefore differ- 
ent representations that are equivalent in the non-oscillating state can differ for the 
oscillating star. It is also evident that the phase lag of the oscillatory perturbation 
to F c behind that of the perturbation to F c q, which controls the damping rate of the 
oscillations, is determined directly by the (arbitrary) value adopted for the constant of 
proportionality between t c and £/w. 

Other approaches address the dynamics of the turbulent perturbations. In the 
mixing-length approach, the fluid is considered to be composed of an ensemble of 
eddies, which somehow attain an existence, grow in the unstable environment and 
then break up, presumably through a ca s cade o f scales. Two extreme approximations 
have been adopted. In one |Goughlll965l ,l 1977al ). the growth is estimated from linear 



theory, thereby ignoring nonlinear interactions, and eventually the nonlinearities are 
posited to act instantaneously to break up the eddy and destroy its correlations be- 
tween velocity and temperature fluctuations. The instants of creation and breakup are 
determined statistically, with probability dist ributions th at are determined solely by 
the background state. In the other extreme ( Unno|[l967l ). the nonlinear interactions 



are presumed to dominate that aspect of the dynamics that determines the turbulent 
fluxes. In a statistically steady state, the eddies are presumed to be steady too, the 
buoyancy driving being balanced by an algebraic representation of nonlinear transfer 
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to other scales. This obviates the need to discuss creation and destruction. Both ap- 
proaches can easily be worked out for a statistically steady convection zone, and yield 
essentially the same formulae for the fluxes. And both can be perturbed for studying 
stellar oscillations, although now the resulting equations for the perturbed fluxes lead 
to very different results. A careful comparison of these two extreme approaches has 
never been undertaken. 

One major deficiency resides in the reduction of the whole turbulent cascade into a 
single length-scale, as it is assumed in mixing-length approaches. As the physical link 
with the cascade is lost, no realistic and reliable predictions for the perturbation of 
the mixing-length can be obtained. Unfortunately, it appears to play a major role in 
mode driving and damping. In all the formulations it is necessary to make assumptions, 
either explicitly or implicitly, about the anisotropy of the convective motion, and how 
it is modulated by the oscillations. In so naive an approach as mixing-length theory, 
this inevitably entails the introduction of new parameters, which, together with the 
parameters in the theory for steady convection, should somehow be calibrated against 
astrophysical observations, laboratory convection or numerical simulation, where pos- 
sible. This is a difficult task, partly because we do not have enough experience with 
modelling and analysing convection and partly because we probably do not have enough 
data in a form in which we can comprehend their implications. The manner in which 
the anisotropy is taken into account affects the formulae for the turbulent heat flux and 
Reynolds stress, and thereby influences the computed damping rates of the oscillations. 

It is important to recognise that the form adopted for the anisotropy of the convec- 
tive motion reappears in the evaluation of the excitation integral of Eq. (|24|l discussed 
in the next section, although it is rare for the formulation of the forcing term 5 in 
the integrand to be consistent with the calculation of the damping rate. The rationale, 
usually unstated, appears to be that because the calculations of both the damping and 
excitation are uncertain, consistency between the two is hardly called for. 

All currently reported stability computations suggest that low-degree solar g modes 
of ord er n > 3 are stable . For these modes, the damping is dominated by radiative 
losses ( Kumar et al.l 19961) and is found to be insensitive to the pulsation-convection 
interaction ( Belkacem et al.ll2009l ). For the stability of g modes with n < 3, however, 
detailed agreement has not been achieved, largely because the balance between nuclear 
driving, radiative damping and energy exchange with the convection is quit e delicate 
( Christensen-Dalsgaard et al.|[l973 ; iBourv et~al1l 19751 ; Ishibahashi et al.lll975h . If grave 
g modes were indeed unstable some nonlinear process must have limited their expo- 
nential growth. 



3.3 Amplitude limitation 



Dziembowskl (| 19821. 1 1983h studied the effect o f three-mode interactions on the g-mode 



amplitudes. In particular. iDziembowskil ( 19831 ) considered the special case of paramet 



ric resonance, in which a parent linearly unstable mode is coupled with two linearly 
damped daughter modes. Once a mode is sufficiently mature, it spawns two resonat- 
ing g-mode daughters who sap energy from their parent at a rate proportional to the 
product of the amplitudes of the three modes. The parent is thereby unable to con- 
tinue to evolve unencumbered, and eventually her growth is stunted. There are various 
ways in which the interactions can proceed, but the most effective for grave low-degree 
parents appears to be a coupling with a pair of similar high-degree twins leading to a 
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state in which the amplitudes of all three modes are steady. The most intricate part of 
the calculation is the evaluation of the interactions. However, the basic overall limit- 
ing process is generic. Therefore we illustrate it with a much simpler, one-dimensional 
acoustic process described in Appendix[X] In that discussion the structure of the back- 
ground state was assumed not to vary with time. Yet to achieve amplitude limitation 
appears to require frequency mismatches Aui of order 10~ e ujQ, which would need to 
be maintained for times of order rj" 1 . That could be several decades or centuries. Can 
such precise resonance be maintained in the face of solar-cycle variation, or any other 
change? And if not, what are the consequences regardi n g g-m ode amplitudes? These 
questions have been addressed by Ijordinson fc Gouehl l|2000l ) , who merely reported 
that the reduction in the amplitude limitation caused by drifting out of resonance is 
offset by the increased probability of actually encountering a resonance, leaving Dziem- 
bowksi's estimate unchanged. 

Dziembowski estimated a likely actual mode amplitude for gi(Z=l) in the pho- 
tosphere to be about 20 cms -1 . That value is considerably greater than the upper 
bounds set by observation. Of course that discrepancy presents no actual serious con- 
flict, because the estimate is probabilistic: Eq. (|5ip does not necessarily preclude a 
precise resonance with daughters of moderate degree whose damping rates are some- 
what smaller than those of the 'most likely' resonance. Moreover, provided solar-cycle 
activity is confined essentially to the convection zone, associated structure variations 
are unlikely to have a large impact on mode resonance, which is determined by condi- 
tions in the region in the radiative interior in which the daughters are trapped. 

But perhaps instead it is more likely that there is no intrinsically unstable g mode. 
All calculations to date have found models of the present Sun to be stable to all g 
modes, although a consistent modern calculation, carried out since the solar structure 
was determined seismologically, appears not to have been reported. So the possibility 
of there being grave self-excited g modes should not yet be wholly abandoned. 

What is it, one naturally might ask, that is so different in more massive stars, such as 
SPB stars or 7 Doradus stars, that enables them to undergo g-mode oscillations at much 
higher amplitudes? The principal difference is that they have only shallow convection 
zones, and therefore daughter g modes extend much higher in the envelope and dissipate 
much more strongly, thus being limited themselves to much lower amplitudes and 
thereby being less able to extract energy from their parent. 

If all g modes were indeed stable today, their amplitudes would be determined 
from a balance between stochastic driving by turbulent convection and damping. We 
therefore discuss in the coming sections the physical mechanisms responsible for solar 
g-mode damping and excitation, and the resulting theoretical amplitude estimates that 
have been published to date. 



3.4 Stochastic excitation by turbulent convection 
3.4.I Amplitudes of oscillations 

If the general solution to the homogeneous equation of free oscillation were known, then 
the solution of Eq. (|22|l for forced oscillations could be written down in terms of The 
general solution for free oscillations of realistic stellar models is not known, but because 
the adiabatic oscillation eigenfunc tions computed in the absence of the int eractions with 
convection form a complete set ( EisenfekH Il969l ; iDvson fc Schutj|l979l ). they can be 
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used as a basis for what is probably a reasonable approximation to the spatial structure 
of the oscillatory motion. What then remains is to characterise the fluctuating driving 

term 5 on the right-hand side of Eq. (|22|l . 

The first serious attempt to address this general problem for the Sun was bv lSteinl 
( 196dll967h . This pioneerin g work addresse d the generation of acoustic waves based 



1ms pioneering work addressed the generation ot acoustic waves based 
on the method developed by iLighthilll (|l952h to derive the acoustic radiation from an 



isolated unstratified free-turbulent region, a region in which no external body force is 
imposed on the fluid. Because acoustic waves can exist in and propagate from the tur- 
bulent region, the local wave-turbulent interaction can be estimated from the turbulent 
field in terms of simple waves. Lighthill showed that, in the relatively straightforward 
situation in which he applied his analysis, conservation conditions required that the 
basic wave emission is quadrupolar. But in a stellar envelope stratified under gravity, 
the situation is more complicated. First, one cannot even define universal multipole 
emission components in general circumstances, although in simple cases, such as an 
i sothermal a tmosphere under constant gravity, a multipolar decomposition is possible 
(|Unndfl96bt ): the stratification produces monopolar (mass fluctuations) and dipolar 
e mission ( force fl uctuat i ons), in addition to quad rupolar (shea r fluctuations), although, 
as lUnnol ( 19661 ). Isteinl ( 19671 ) and subsequently lOsakil (|l99Cj ) pointed out, there is a 



tendency for the monopolar and dipolar components to cancel. 

A direct generalisation of Lighthill's method is not possible, because gravity waves 
cannot propagate through a convectively unstable medium. Indeed, it is probably not 
most prudent even to work in terms of simple emitted waves, because one would still 
have to impose resonance conditions on them in order to represent the normal modes, 
which, as Stein appreciated, involves some subtle argument. Instead, it is more straight- 
forward to solve Eq. (|22p as a nonsingular perturbation, regarding 5 as a forcing term 
on a spatially unperturbed oscillation eigenfunction £ of oscillation having angular fre- 
quency u) and (positive) damping rate 77. Basically, one represents the forced motion as 
A(t)£(r) exp(— itot — rji), where A varies slowly compared with exp(— \ujt). Substituting 
this into the spatio-temporal wave operator from which C was derived, and retaining 
only terms proportional to A and its first derivative, yields a differential equation for 
A whose solution can be written down explicitly as an integral depending on £ and g 
(iGoldreich fc Keelevlll977bl ; lBalmforthlli992al : iMusielak et al]|l994l : ISamadi & : Goupill 
l200ll ; IChaplin et al.ll2005l ). For statistically stationary turbulence exciting modes hav- 
ing characteristic wavelengths much greater than the typical scale A; -1 of all convective 
eddies and with r\ « to, the integral expression for the mean-square amplitude A re- 
duces to 

1 r°° 

A 2 ~ / G(r) <T VT cos {lot) At (24) 



ir/cj 2 ! 2 

( Chaplin et al"ll2005h . where G(r) = Q(t')Q*(t") is a global correlation function with 



argument t = t — t , the overbar denoting ensemble (or temporal) average, and 

Q(t) = I C(r)-5(r,t)dV, (25) 



' V 

the integral being over the volume V of the star. The asterisk denotes complex conju- 



gate, and 1 = J p^.^dV is the inertia of the mode. The ensemble average Q(t')Q*(t") 
can itself be written as a double integral over the volume of the star. Q(t) depends on 
a two-point correlation function between (r' , t') and (r" ,t") which can be represented 
by its (local) Fourier transform x{ u > r ) with respect to r and r — r referred to the 
mean position r = (r + r")/2. 
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The principal uncertainties in the estimation of A 2 lie in the evaluation of the 
damping rate 77 and the correlation function G. The mean-squared surface velocity is 
computed, for each mode, as 

^j = A^[v?(h)+vl(h)] (26) 

where h is the height in the atmosphere at whic h the complex oscil l ation velocity 
(vi,V h) = iwo(£r,£h) i s observed (for details see IHoudek et all 1 19951 : IBaudin et al.l 
120051 ), 



3.4-2 Application to g-mode amplitudes 



There have been several attempts to estimate the stochastic excitation of g modes, 
all of which stem from the early formulation by Goldreich & Keeley (1977a,b) of the 
mechanism of t h e exc itation of global p modes. This work was initially extended by 
iGoldreich et~aT] (|l994T ) and then by others. The analyses differ in detail, as do the 
results. The differences in the assumptions adopted in the analyses occur in several 
places, and it is not possible to determine from the published work just how they 
influence the conclusions. Therefore we shall only summarise very briefly the principles 
behind the work. 

The intention of Goldreich and his col leagues was to impro ve the earlier work on the 
excitation of solar p modes. Subseauentlv lKumar et aD ( 19961 ) addressed the issue of g- 
mode observability. The first step in the endeavour was to derive an explicit form of the 
linearized momentum equation governing the oscillations in order to obtain an explicit 
expression for the forcing term in Eq. (|22[1 . IGoldreich et al. ( 1994 ) approached that 
task by implicitly assuming that the convection also satisfies the linearized momentum 
equation. The outcome is 



V • (pvv — pvv) 



(27) 



where p is pressure, although IGoldreich et al.l (| 19941 ) did not explicitly separate the 
fluctuations from th e mean s. This result generalises a similar expression presented by 
IGoldreich fc Kumarl (|l990l > for an isentropically stratified plane-parallel atmosphere of 
perfect gas. The function 5 was entered into the integrals in Eq. (|25fl and the terms 

estimated from mixing-length theory. 

In deriving the inhomogeneous wave ecmation lGoldreich et al.l ( 1994j) omitted terms 
containing the Eulerian density perturbation (even when coupled with the gravitational 
acceleration, a term that must not be neglected even in the Boussinesq approximation) 
and the inertia terms of the convection. Furthermore, the authors assumed that the 
effect of the oscillations is to destroy the balance of the convective momentum trans- 
fer in such a way that the linearized momentum equation is satisfied in the absence 
of Eulerian entropy fluctuations, leaving the unbalanced convective entropy fluctua- 
tions and nonlinear inertia terms to drive the oscillations. This assumption caused the 
entropy fluctuations to be an order of ma gnitude more effective than the fluctuating 
Reynolds stress in driving the oscillations. IGoldreich et al.l ( 19941 ) ignored products of 
oscillation variables and convective fluctuations, and the effect of the oscillations on the 
convection. Consequently an adiabatic linear wave equation (with homogeneous terms 
untrammelled by convection) is supplied just with an inhomogeneous forcing term to 
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account for all the effect of convection. The anisotropy of the convective motion was 
acknowledged merely by multiplying the entire inhomogeneous term S by a single con- 
stant scaling factor, which was determined by calibrating the energy supply rate to the 

oscillations against observation. 

The theory wa s reformulate d bv lSamadi fc Goupill d200lh. who combined aspects 
of the analyses by Isteinl 1 1967t ) and IGoldreich fe Keelevl i|l977bl ). Perturbations due 
to convection and global oscillations were separated in the equation of motion, and 
approximated according to their dominating dynamics, a dopting a Boussinesq-like ap- 
proximation for the convection. ISamadi fc Goupill ( 200ll) have shown th a t the linear 
term due to entropy fluctuation [Eq. (|27[l ] introduced by IGoldreich et ah I (| 1994 ) gives 
no significant contribution to the driving but the advection of Eulerian entropy fluctu- 
ations by turbulent velocity. The analysis represents what is probably the most careful 
recent attempt to generate a consistent forced wave equation. Unfortunately, however, 
there is subsequently an error in the treatment of the frequency correlati ons in the ex- 
pressi on for the driving by the fluctuating Reynolds stress, corrected bv lSamadi et alj 
(|2Q05h . 

The pr edictions for the max imum amplitude for £ Hydrae made bv lHoudek fc Gough 



20021) and lSamadi et all (|2007T l are in good agreement with the observations by; 



Frandsen et al.l 



20021 ) . while their predicted maximum amplitude for Procyon are overestimated by 



a factor 2 to 4. This discrepancy for Procyon is serious and must be understood; the 
structure of Procyon is rather different from other stars that have been modelled, in- 
dicating that the theory, anchored by calibration against solar p modes, cannot not 
be extrapolated reliably to very different stars. By the same token the reliability of 
extrapolating from p modes to g mod es even in the Sun m ust be exercised with due 
caution, and may not be as reliable as lKumar et al. (1996) s uspect. 

One of the controversial conclusions of the discussion bv IGoldreich et all ( 19941 ) is 
that driving from the first term in Eq. (|27|l for 5 is an order of magnitude greater 
than driving from the fluctuations in the Reynolds stress, at least fo r acoustic emis- 
sion. This contradicts the earlier findings of lGoldreich fc Kumarl 1 1990l ) who concluded 
that, in a convection zone that appears to the oscillations (but not to the convection) to 
be adiabatically stratified, the contribution from the two terms are comparable. The 
latter conclusion ha d been drawn from the realisation that the emission from what 
IGoldreich fc Kumarl (|l990h termed the monopole and dipole sources (con tributions 
which in a stratified envelope are not unambiguously defined, cf. lUnnol 1 19661 ) from the 
buoyancy (entropy fluctuation) terms largely cancel as a result of the manner in which 
the eddy motion is correl ated w i th the m. Such cancellation r esults in predominantly 
quadrupolar emission, as lOsakil (|l990l . see also iHoudekl l2006h also pointed out, and 
is consis tent with Stein's (1967) e a rlier finding s . It is interesting to note that simula- 
tions bv IStein fc Nordlundl d200lh; IStein et ail jjooj ), and the model calculations by 
iBelkacem et ail <|2006bl ) and lSamadi et ail (|2007m indicate that emission is dominated 
by the fluctuating Reynolds stress and entropy fluctuations as had been implied by 
Balmforth's (1992b) discussion. 

The treatment of the turbulent velocity spectrum in evaluating th e correlation inte- 
gral of Eq. (|24|l deserves some comment. In this context. Isteinl ( 19671 ) first adopted the 
quasi-normal approximation to relate the fourth-order correlation to a sum of three 
independent products of second-order correlations, as have othe rs who have explic - 
itly considered the velocity correlation after him, notwithstanding iKraichnan'sl (|l957l ) 
warning of the danger of so doing. This decomposition would be correct if the probabil- 
ity distribution function of turbulent quantities were Gaussian, which is not what was 
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subsequently assumed fsee lBelkacem et al.ll2006al : iKupka fc Ro binson 20071). Note that 
a decomposition of the four th-order momen t that t akes the effect of plumes into ac- 
count has been pr oposed by Belkacem et alj ( 2006al ) and permits to better reproduce 



the observations (jBelkacem et al.1 l2006bl ). The dominant term (a product of single- 
point correlations) is balanced by the Reynolds stress that appears as the last term 
in Eq. (|27[) . leaving the forcing t erm t o dep end only on two-point correlations. The 



next approximation adopted by Isteinl ( 19671 ) is to write the energy spectrum as a 
product of a function of local wavenumber and a function of frequency, relating the 
three-dimensional velocity correlation to the energ y spectrum by th e formula valid for 
incompressible, homogeneous isotropic turbulence (|Batchelorlll953h . This result is ap- 
proximately valid for isotropic turbulence in the Boussinesq approximation - almost 
universally adopted, either explicitly or implicitly, in mixing-length theory - and which 
can be g eneralised to the ax isymmetric case expected of convection in a nonrotat- 
ing star ( Chaplin et alj 2005), although the generalisation to anisotropic turbulence 
has been adopted fe. g. iBelkacem et alj 20091). Similar appro x imati ons h ave been used 



p-mo de studies by iBalmforthf (|l992al ). ISamadi fc Goupi] |200ll ) and ISamadi et alj 
(|2003bh . 



3.4-3 The energy- equipartition principle 

A simple way of crudely estimating the mode amplitudes, without recourse to an 
explicit model for the excitation and damping of oscillat i on mo des, is to adopt the 
equipartition ansatz discussed first bv lGoldreich fc Keelevl (|l977bl ). It is obtained from 
a rough estimate of the integral in Eq. (|24|l . which is proportional to the ratio of the 
rate of forcing to the energy lost by dissipation (cf iBatchelorl fl953l ) . The energy sup- 
ply rate was estimated from a mixing-length-like description of turbulence assuming 
that the motions of all the convective eddies are independent of one another, and 
that, as usual, the spatial scale of the eddies is the smallest of all scales, being less 
than both the scale height of the background state and the inverse wavenumber of the 
oscillation eigenfunction. The temporal spectrum of the turbulence was taken to be 
Gaussian, with variance proportional to the square of the characteristic timescale r c of 
the energy-bearing eddies (i.e. the largest eddies, with spatial scale equal to the mixing 
length I), and with an autocorrelation that is significant only o ver the spatial scale 
of a s ingle eddy. The intrinsic damping rate n was estimated bv lGoldreich fc Keelev 



1977 ) from a radial-pulsation calculation in which they presumed, as had lCox et al 



19661 ) previously, that the phase of the perturbed heat flux F{. is determined by the 
equation dF^/dt = —F c /t c , and that the Reynolds stress, whose influence dominates 
the energy loss, can be represented by a time-independent scalar turbulent viscosity 
v% whose value is either v^t c , where v c is a characteristic velocity of an energy-bearing 
eddy, provided lut c < 1, or v\t\, where v\ is the velocity of an eddy in the turbulent 
cascade (assumed to satisfy Kolmogorov scaling) that resonates with the pulsation: 
ojt\ = 1. This discont i nuous change in behaviour had been adopted previously by 
IGoldreich fc Nicholsonl (|l977t ). and takes some account of the fact that eddies with 
timescales much greater than lo~ 1 contribute very little to the dissipation of th e pulsa- 
tional motion; more soph isticated mixing-length descriptions of convection (e.g. iGoughl 

Il965l . Il977al ; IUnnolll967l ) predict a c ontinuous transition. 

From simple scaling arguments IGoldreich fc Keelevl ( 1977bl ) concluded that the 
turbulent fluctuating momentum fluxes with timescales co mparable w ith the pulsa- 
tion period dominate the driving term Q [Eq. (|25|l ], as had Isteinl ( 19671 ) before them, 
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and, because the dissipation is also dominated by the same fluxes, via ut, there is a 
cancellation of a squared turbulent velocity in the power-to-dissipation ratio, leading 
to an equipartition between the energy 7? sc in a mode of oscillation and the energy 
p£ 3 Vc/2 in a single resonating eddy satisfying lut c = 1, if such an eddy exists, or the 
most energetic resonating eddy in the turbulent cascade if no such energy-bearing eddy 
exists. In the case where one does exist, this balance may be written 

E osc = A Pt l 3 = E c , (28) 

where pt is the (r, r) component of the Reynolds stress (also known as the turbulent 
pressure), and A is a factor of order unity whose value depends on uncertain properties 
of the turbulent flow, and is therefore itself uncert ain. 

Th is equipartition principle was adopted by IChristensen-Dalseaard fc Frandser] 
(1983) for estimating p-mode amplitudes in main-sequence stars ; and was applied for es- 



timat ing g mode amplitudes in the Sun by Gough (1985; see also lBerthomieu fc Provostl 
1 19901 ;). In these calculations it was necessary to estimate the convective timescale 
r e , which in the Su n has a minimum in the superadiabatic boundary layer (see e.g. 
IChaplin et ailbood ). Hence, for a mode with a period 77 that is longer than min(r c ) 
there are at least two radii in the convection zone where r c = 77, because formally r c 
tends to infinity towards the boundaries of the convection zone. The kinetic energy in 
a mode is then balanced against the sum of the values of E c at these points. 

Because the value of A is uncertain, iGoughl (|l985l ) calibrated the calculation by 
equating the maximum value of the amplitudes of the 5-minute p-mode oscillations 
with the value observed. He then obtained a maximum velocity amplitude of about 
0.5 mm s" 1 for the gravest quadrupole g mode and lmms -1 for the gravest dipole g 
mode, and about 3 mms -1 for gravest of low-degree p modes. The result is very nearly 
a function of frequency alone; it is depicted in Fig. 1151 

3.4.4 The results of Kumar, Quataert and Bahcall 



Kumar et ail J 19961) . 



jThomson et all 



iGoldreich et al 



mulation. 



motivated by a claim of g-mode detection in the solar wind 
1995), carried out computations using the formalism developed by 



( 1994 ) origin ally to estimate the amplitudes of p modes. In this for- 



Kumar et al] ( 19961 ) assumed a simplified description of turbulence in which 



the velocities, length scales and timescales cascade consistently according to the Kol- 
mogorov spectrum. They assumed that the modes are dri ven solely by the fluctu ating 
Reynolds stress. Their calculation was identical to that of IGoldreich et aL I jl994h . 



Of particular interest is t he way in which the eddies and the standing waves are 
temporally correlated. In the Goldreich fc Keelev (1 1977b ) approach , from which was 



derived the procedures of iBalmforthl ( 1992a! ) ; IGoldreich etall (|l994T ). it was assumed 



that the temporal correlation between eddies is Gaussian. As we shall see in the fol- 
lowing sections, variations in the manner in which that function is chosen can lead to 
very different estimations of g-mode amplitudes. 

Concerning the damping rates, both turb ulent and radiative c o ntribu tions to the 
damping rates were included as derived by IGoldreich fc Kumarl |l99ll ). Turbulent 



damping was approximated purely as a momentum-diffusion process, with a time- 
dependent scalar diffusion coefficient derived from mixing-length estimates of the mo- 
tion of eddies with lenghtscales smaller than the local wavelength of the oscillation. 
Temporal modulation of the Reynolds stresses, considered to be as i mportant as mo- 
mentum diffusion, if not more so, at least for grave p modes (e.g., iBaker fc Goughl 
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Fig. 14 Theoretical damping rat e r\ of g modes of low degree plotted versus cyclic frequency 
v according to the predictions by Bclkaccm ct al. ( 2 009f l . 



Il979l ; lBalmforthlll992bh . was ignored, as was the modulation of the energy generation 
in the core. For higher-order p modes of low degree, the diffusive contribution was 
estimated by iGoldreich fc Kumarl ( 199ll ) to be comparable to the radiative losses, bu t 
for low-order p and g modes it is smaller (|Goldreich fc Kumar|[l99ll ; lBalmforthlll992bl ). 

IGoldreich fc Kumarl i|l99lf) found mode lifetimes of about 10 6 y, comparable with 
previous estimates. The resulting surface mode amplitude, based on Eq. (|26[) , was found 
to be greatest, at about lmms -1 , for I = 1 modes near v — uj/2ir 



with v < 100 /iHz were all found to have amplitudes less than 10 mms 



200 nRz; modes 
1 



3.4-5 The results of Belkacem et al. 



iBelkacem et all (|2009h investigate d the particular case o f the amplitude of asymptotic 
g modes. The formalism used by IBelkacem et aL ( 200S ) to co mpute excitation rates 
of non-radial modes was developed by IBelkacem et al (|2008h who extends to non- 
radial modes the work of ISamadi &i Goupil (2001). The turbulent-mean stresses were 



neglected since the excitation takes place in the innermost convective layers where the 
Mach number is very small. The driving is shown to be dominated by the Reynolds- 
stress fluctuations. In this formalism, the turbulent kinetic energy spectrum as well 
as the eddy-time cor relation function are derived from 3D numerical simulation (see 
IBelkacem et al.ll2009l . for details). The Lorentzian function (xk) ls found to better 
reproduce the eddy-time correlation function from the 3-D numerical simulation than 
a Gaussian function in the frequency range v 6 [20 ^iHz; 110 ^iHz]. In addition, the 
eddy-time correlation function is poorly represented by a Gaussian function, which 
underestimates \k by many order of magnitudes. 

Da mping rates are com puted with a fully non-radial non-adiabatic pulsation code 
MAD jDupret et al.ll2006d lbllal) , including a non-local tr eatment of convection based on 
lUnno (|l967T ) formalism extended to non-radial modes bv lGabriell |l996l ) ; iGrigah cene et al.l 
(2005). It takes into account the role played by the variations of the convective flux, the 
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Fig. 15 Estimated amplitudes of stochastically excited g modes of low degree plotted against 
cyclic frequency v. The estimates are rms surfac e values for s inglet modes (si ngle values of 
re, I an d m), joined b y lines: triple-dqt-dash ed for iGoueh (1985J), dot-dashed for lKumar et all 
1 119961) and dashed for Bclkacem et al. N2009D . The continuous line is an estimated from 10 years 
GOLF data l|Garci'a et al.ll2007l V For prop er comparison, the eff ect of the spatial instrumental 
filter must be included, i.e. visibilities (Sec Belkacem et al. 2009) (Sec Section on Obscrvables). 



turbu lent pressure and the dissipation rate of turbulent kinetic energy. iBelkacem et alj 
(2009) have found that for high frequency g modes (y > 110 /iHz), the work integrals 



and thus the damping rates are sensitive to the convection/pulsation interactions be- 
cause the role of the surface layers in the work integrals becomes important. In contrast, 
for low-frequency g modes [y < 110 /^Hz), the work integral and then the damping rates 
are found to be insensitive to the convection/pulsation interactions as well as the non- 
local parameters. Eventually, the damping rates dominated by radiative losses behaves 
„o ,,-3 



Ta king visibility factors as well as the limb-darkening into account. IBelkacem et al.l 
( 20091 ) inally found that the maximum of apparent surface velocities of asymptotic g- 



modes is w 3 mms 1 for £ = 1 at v ~ 60 /iHz and i = 2 at v ~ 100 ^iHz. Those results 
then put the theoretical g-mode amplitudes near the GOLF observational threshold. 

The principal differences between the calculations lBelkacem et al] ( 20091 ) and lKumar et all 



( 19961 ) stems from the turbulence spectra that were adopted. It would create enormous 



differences in g-mode amplitudes. They are produced by convection-oscillation interac- 
tions deep in the convection zone, which are off resonance, and whose intensity depends 
crucially on th e assumed fr eque ncy dependence of the turbu lent spectrum in the high- 
frequency tail. ISteinl (| 19671) and lGoldreich fc Keelevl (Il977bl) had favoured a Gaussian 



distribution. A Gaussian was preferred also by Chaplin et al.l (|2005t ). for although the 
energy-bearing eddies of the convection theory they used (that of Gough 1965, 1977a) 
has a Lorentzian time-dependence. We discuss the disparity in Section \3. 51 
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3.4-6 High-order waves generated at the base of the convection zone 

A series of papers using a wave-mechanical consistent numerical approach with a sim- 
plified two-dimensional solar model attempted to quantify the excitation of gravity 
waves in the solar in terior as well the transmission of gravity waves through the so- 



lar convection zone (Press 



1981 



Zahn[ 1991 Andreas sen et all 1 19921 : lAndersen et al.l 



19931 : lAndersen! Il994l . , 1 1990 : iTalon fc Charbonne]||2003l ). The results indicate that the 



efficiency of generating gravity waves in the interior by convective overshoot is slightly 
less than 0.1% for waves with a horizontal extent equivalent to 1=6-18. The overshoot 
layer where the excitation occurs is very shallow, 0.02-0.06 pressure scale heights. The 
energy transmission through the convection zone of a gravity wave equivalent to an 
1=6 mode is found to be 3-8 10~ 4 in the 50-250 /iHz. Combining these results indicate 
that the convective overshoot into the interior provides enough energy to sustain an 
amplitude of 1-5 mm s~ at the surface of a single gravity mode. If the energy is dis- 
tributed into more modes the amplitudes are reduced by the square root of the number 
of modes. 

These results should not be taken as serious physical predictions of amplitudes of 
solar g-modes, but rather as a demonstration that there is enough energy available in 
the convective overshoot to maintain a single g-mode at an observable level at the solar 
surface. Assuming that a thousand modes are excited and scaling the amplitudes to an 
1=1 mode implies an amplitude of about 0.02-0.08 cm s^ 1 in th e 50-250 uHz. As seen 
from Fig. 1151 this is approximately the same values provided bv lGoughl (f 19851 ). 



3.4-7 Prediction of g-mode amplitudes 



Amplitudes o f singl et g modes (single values of n, I an d m) predicted bv lGoughl ( 19851 ) 



Kumar et ail (119961 ) and iBelkacem et all (|2008i . 120091 ) are depicted in Fig. Q2] Shown 



are the root- mean-square total velocity amplitudes in the photospher e; instrumental 
filtering (e.g. lDziembowskill 19771 ; [christensen-Dalsgaard fc Goughlll982l ) must be taken 
into account when converting these values into observables. 



3.5 Discussion 



The p rincipal difference between the calculations of lKumar et alJ ( 19961 ) and lBelkacem et al.l 
(|2009h stems from the turbulence spectra that were adopted and, in particular, the way 
the two lead to enormous differences in g-mode amplitudes. Those differences are pro- 
duced by convection-oscillation interactions deep in the convection zone, which are off 
resonance, and whose strength depends critically on the assumed frequency dependence 

of the turbulent spect rum in the h igh-f requency taiffj 

For solar p modes. Isteinl ( 19671 ) and lGoldreich fc Keelevl ( 1977b) had used a Gaus - 
sian time-correlation function. A Gaussian was preferred also by Chaplin et al.l ( 20051 ). 
for although the energy-bearing eddies of the convect ion theory they used (that of 
Gough 1965, 1977) has a Lorentzian time-dependence, IChaplin et al. I |2005l ) reported 



5 In this section, we denote by time-correlation function the Fourier transform Xk of the local 
correlation func tion entering the evaluation of the global co rrelation function G(t) appearing 
in Eq. 1 1241 (sec Samadi & Goupil 2001; Chaplin ct al. 2005, for details). We discuss Gaussian 
and Lorentzian forms, both of which have characteristic widths that depend on the lengthscales 
and timescales of the convective eddies, and therefore vary with position r in the star. 
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that th ere is experimental evidence for more rapid decline in the tail. However, Samadi et alj 



2003al had pointed out that the three-dim ensional simulations by IStein fc Nordlun 



200 ll w ere more nearly Lorentzian. Indeed. ISamadi et al] ( 2003aT ) and 



Bclkacem et al.l 



( 2006bl ) reported that the use of a Lorentzian function repr oduces the observati ons 



for solar p modes better. A similar conclus ion was reached by Samadi et al] (|2008f l in 
the case of the star a Cen A. In contrast, IChaplin et ail (fcooa i" mention that a pure 
Lorentzian results in an over-estimation for the low-frequency modes. They explained 
that the over-estimation of the mode excitation rates at low frequency results from 
the Lorentzian function decreasing too slowly with depth compared to a Gaussian. 
Consequently, a substantial fraction of the excitation rate of the low-frequency modes 
arises from large eddies situated deep in the Sun. A Gaussian time-correlation function 
gives much less weight to the large off-resonance eddies. Which is more nearly correct 
remains an open issue. 

iBelkacem et al] (|2009h favoured a Lorentzian modelling for solar g modes, on the 
groun d that it reproduces the 3D numerical convection simulations by iMiesch et al] 
( 2008) more closely. More specifically, the best fit was found using a sum of Lorentzian 
and Gaussian functions. The Gaussian shape is found by the authors to reproduce 
the very low- frequencies while the Lorentzian the highest frequencies, which are of 
particular interest. There is yet no definite issue, but what is sure is that it has a 
crucia l impact on predicted g-mode amplitudes. As demonstrated bv IBelkacem et all 
( 20091 ). Gaussian and Lorentzian time-correlation functions lead to differences in mode 
surface velocities of several order of magnitude. 

It is worthwhile mentioning that experimental studies of turbulent convection can 
provide some c l ues to wards the resolution of this issue. In particular, the recent work of 
iMordant et al] ( 20041 ) demonstrates that the two-point second-order time- correlation 
relevant to our discussion is reasonably well reproduced by an exponential decrease 
with separation t i me in the inertial range (see t heir Fig. 10) w hich is in agreement 
with iKolmogorovl ( 194ll ) phenomenology (see also ISawfordlll99~lh . It corresponds to a 
Lorentzian power spectrum with respect to frequency, which lends support to the use 
of a Lorentzian in the Sun, even though the Prandtl number in the experiment (6.8) 
is very far from solar. At frequencies beyond th e inertial range dissi pation causes the 
spectrum to decline more rapi dly (see Fig 13 of IMordant et all I2004T ) . It is interesting 
to record that simulations by iGeorgobiani et aL |2006h ind icate that the spec trum 
is actually not separable in frequency and wavenumber, as iKraichnan'sl ( 1957t ) had 
already pointed out, and suggested an empirical improvement. 

Another import ant matter for solar g modes is the way in which the damping 
rates are computed. IBelkacem et all ( 20091 ) did not calibrate their predictions against 
observations of p modes, because they considered that because the g modes of higher- 
frequency (y > 110 Hz), whose damping rates depend on the uncertain modelling of 
the interaction with convection, could not be predicted reliably, a meaningful connec- 
tion with the p-mode regime of current observation could not be made. Instead they 
restricted attention to low-frequency modes whose damping rates r\ are dominated by 
radiative losses, and scale with u~ 3 , a behaviou r which follows from straightforward 
asymptotic argument fe.g. IBelkacem et alj |2009T). and which is evident in Fig. 1141 But 
for grave g modes, the situation is less clear. It should perhaps be pointed out in 
passing that the e mechanism is significant for the grave g modes, and the influence 
of the dynamically induced thermodynamic perturbations must consistently be taken 
into account. A predictive description of the interaction between convection and oscil- 
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lations when both are strongly coupled is mandatory to obtain reliable mode amplitude 
estimates in this frequency domain. 



3.6 Some concluding remarks 

Concerning the damping rates, the failure to detect g modes f avours that they are 
linea rly st able. Theoretical amplitud e estimates investigated by iDziembowskil f|l983l . 
1982) and Ijordinson fc Goughl |2000l ) are well above the current detection threshold, 



although it must be appreciated that the results are somewhat uncertain. 

In the high-order asymptotic regime, radiative diffusion is thought to be the domi- 
nant contribution to g-mode damping. For low-order g modes the situation is less clear 
since those modes are sensitive to the interaction with convection, as are the p modes; 
there is no co nsensus about the dominant contribution to the mode line- width (e.g., 
iHoudeyEooj ). This issue is very important and prevents an unambiguous theoretical 
determination of low-order g-mode amplitudes. 

For the drivin g mechanism, turbulent convection is thought to be r esponsible for 
g-mode excitation (|Gougblll98 O nKumar et al.|[l996l ; lBelkacem et al-lboolil ) . and is prob- 
ably dominated by the contribution of the Reynolds stresses. Indeed, the so-called en- 
tropy contribution to the driving is found in the current theoretical computation to be 
negligible for both low- and high-order g modes since they are not sensitive to the very 
upper layers where entropy fluctuations are important. 

Quantitative estimates of mode amplitudes differ from each other by orders of 
magni tude, depending pri n cipally on the way th e turbulent eddies are time-correlated. 



While iKumar et ail (Il996l);[chaplin et al l d200otl favou rs a Gaussian p rofile to describe 
this fu nction, ISamadi fc Goupill l|200lh ; ISamadi et all (|2003al . 120071 ): iBelkacem et all 



( 20091 ) used a Lorentzian profile. It should be clear that the eddy correlation function 



cannot be analytically derived, and therefore is adhoc taken either from 3D simulations 
or turbulence experiment; as such it is a scale parameter. In fact the truth is probably 
somewhere in between, the challenge is now to determine where exactly. The detection 
of g modes would provide a definite answer to this burning question. 



4 Detection techniques 

In this section, we present different detection techniques that have been applied in 
searches for g modes. We discuss the conceptual design of each technique using the 
metaphor of bricks. We consider each detection technique to be built or comprised of 
several bricks or categories. An appropriate combination of various bricks constitutes 
a given detection technique. The categories we consider are as follows: 

— Observables 

— Spatial and temporal filters 

— Spectrum estimators 

— Statistical testing 

— Patterns 

— Data combinations 

The choice of the solar observable (or perturbation) to be observed plays a key role 
in determining the likely success, or otherwise, of a given technique (i.e., it can affect 
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sources of noise, observed mode amplitudes etc.). Since we analyse time series, the 
temporal filter of these series (e.g., low-pass, high-pass, band-pass) also merits careful 
attention . Observations of the Sun may be made by imaging the surface onto several or 
many detector pixels, or by observing the "Sun as a star" . When resolved observations 
are made, it is advisable to apply appropriate spatial masks to extract the signatures 
of the modes of interest. 

The detection of potential frequencies of interest requires that we obtain estimates 
of the frequency spectrum, using Fourier transform or other techniques. Then, a thor- 
ough statistical assessment of what is observed plays a key part in deciding whether 
or not something of interest has been detected. If the identification of individual fre- 
quencies fails, one can then bring more a priori information to bear in the detection 
procedures, for example by searching for signatures of patterns expected from theoret- 
ical modelling of the g modes (e.g., near constant spacings in period of modes in the 
asymptotic regime, rotational frequency splittings etc.). Use of a priori knowledge in 
this way suggests a Bayesian approach to the analysis. Finally, if the detection is not 
successful with a single instrument, one can bring data from different instruments to 
bear, for example by searching for coincidences in the different datasets. 



4.1 Observables 

Perturbed parameters such as temperature and displacement can be deduced from 
observables. The displacement can be either physically observed or inferred from the 
solar radial velocities. 

The solar oscillat ions were first detected in radial velocity observations made by 
ILeighton et~aTl l|l962h . This method is now widely used for observing solar p modes. So- 



lar ra dial velocities can either be measured from resolved images of the Sun (|Unno et al.l 
1989, and references therein), or from full-disc integrated Sun-as-a-star observations. 
The instrument aboard SOHO all rely on measuring the intensity distribution in an 
absorption line for deri ving the solar radia l velocity either locally such as with the 
SOI/MP^ instrume nt (|Scherrer et al.lll995l ) or globally with the GOLff] instrument 



( Gabriel et aHll995l ). This set of instrument is also co mplemented by grou nd-based 
networks such as the high-resolution GO NG3 instrument dHarvev et a"flll996l ) , and the 
full-disk integrated BiSON^! instrument ijChaplin et"al]|l996t ). The lengths of the us- 



able time series are about 14 years for MDI and GOLF, 15 years for GONG and nearly 
30 years for BiSON. 

Temperature perturbations induce intensity fluctuations that can be detected. In- 
tensity fluctuations may also be induced by c hanges in the size an d shape of the so- 
lar surface produced by the oscillation modes (|Toutain et al.lll999h . Intensity fluctua- 
tions were first observed from _space by IWoodard fc Hudson! (119831) in solar irradiance 
data collected by the ACRIlvfj instrument |Willsonlll979l ). Simi lar observation s have 
been made with sunphotometers from stratospheric balloons (e.g. lFrohlich1ll984l ) which 



confirmed the ACRIM results. Ground-based observations with the SLOT] 11 ! instru- 
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merit produced rather marginal p-mode detections ( Jimenez et al.lE987l ). while rather 
higher signal-to-n oise data were colle c ted by the LoJ 12 l instrument, observing from 
El Teide, Tenerife lAppourchaux et al, | |l997l , ll995l ). Following the success of ACRIM, 
the IPHIe(3 instrument onboard the PHOBOS mission lead to the firs t long-duration 
intensity measurements of the solar p modes ( Toutain fc FrohlictJ 19921 ) . The ACRIM, 
SLOT and IPHIR instruments lead to the design of the VIRG(£3 

instrument for the 

SOHO mission, which was composed o f two radiometers (D IARAD and PM06), three 
Sun PhotoMeters (SPM) and the LOI |Frohlich et al.lll995l ). The lengths of the usable 
time series are about 14 years for VIRGO. 

Intensity fluctuations produce a much lower signal-to-noise ratio in the p modes 
than do solar radial velocity observations. The signal-to- noise ratio for the former is at 
most 30, while the latter is at most 300; exa mples of the various signal-to-noise ratios 
for the different observables can be found in iToutain et all (|l997T ). 

The detection of the physical displacement is more difficult. For instance a putative 
g-mode amplitude of 1 mm s^ 1 at a frequency of 100 /xHz would induce a peak-to-peak 
variation of 4.3 /^arcsec as seen from the Earth, whereas for a p-mode amplitude of 
30cms _1 (at a frequency of 3000 /iHz) would induce a variation of 43 /^arcsec. In the 
1970's, the SCLERA^ in strum e nt was used in an attempt to detect p modes and g 
modes alike |01eson et al J 1 1974 lHilj||l985l . 1 19921 ) . The detec tion of displacemen t was 
performed by observing the change in diameter of the Sun ( Brown et al. I [19781 ). and 
signals with typical amplitude of the order of 5 marcsec were observed. These ampli- 
tudes were about a factor of 100 larger than the expected amplitudes for p modes, 
based o n high signal-to-noise detections of the p modes made in radial velocity obser- 
vations ( Claverie et al.lfl98ll : iFossat et al.lll98ll ). This discrepancy cast some doubt on 
the identification of the detected SCLERA peaks as p modes. 

The PICARD mission, due for launch in late 2009 (jPame et all 1 19991 ). will be 
sensitive not only to the displacement, but also to intensity fluctuations induc ed by 
temperature perturbatio ns, like those detect ed bv lAppourchaux fc Toutainl ( 19981 ) with 
the VIRGO/LOI or by iToner et all l|l999l ) with SOI/MDI. PICARD will have the 
capability of making observations of displacement and intensity fluctuations at the 
limb of the Sun, where the perturbations appear to h ave higher amp l itude than those 
observed in the disc. This amplification, theorised bv I Toutain et all ( 19991 ). is due to 
the fact that at the solar limb the atmosphere becomes more transparent, providing a 
larger contribution from the other hemisphere. 

Another potential observable is the perturbation of the gravitational field caused 
by the g modes producing tida l perturbations of the Newtonian fields , and generating 
detec table gravitational waves (jChristensen-Dalsgaard fc Goughlll98b1 ; lGiamperi et all 
2000). The perturbation is non-zero only for modes with / > 2; for I = 1 it is ex- 
actly zero. The advantage is that the solar noise caused by convection would have 
a negligible impact on the pe rturbation of the gravitational field. Unfortunately, as 
IChristensen-Dalsgaardl (|2002bl ) explicitly mentions, signal from galactic binary stars 
could be the limiting factor for g-mode detection and not the Sun itself. 
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4.2 Spatial and Temporal filters 
4-2.1 Temporal filters 

The analysis of discretized signals dates back to the start of the telephone and tele- 
graph ic system, when Harry Nyquist devised the so-called Nyquist frequency (|Nyquistl 
1924). and when Claude Shannon enunciated the Shannon theorem concerning the 
sampling of band-limited signals ( Shannonlll949h . Both worked at Bell Telephone Lab- 



oratories. 

Although, we do not always follow the Shannon-Nyquist prescription, one should 
not forget the limitation inherent in the sampling of a band limited signal. The main 
difficulties encountered in astrophysical applications are as follows: 

— the signal is not band limited 

— the signal is undersampled 

— the signal is aliased 

These three factors are not unrelated to one another. 

In our case, the observations of intensity fl uctuations and of solar radial velocities 
do n ot present a clear high-frequency cut-off ( Kosovichev et all 1 19971 ; iFrohlich et al.l 



19971 ). The freq uency drop-off of the solar backgro und noise at high frequencies is pro 



portional to v z (lHarvevlll985 ; Aigrain et al.ll2004]. a nd references therein), or possibly 
even a higher power law |Appourchaux et alj|2002h . Since the astrophysical signal is 
not band limited, the highest frequencies will not be properly sampled, resulting in an 
undersampled signal. This leads to high-frequency signal leaking back into the spec- 
trum, below the Nyquist frequency, i.e., aliasing of the signal. The imaging effect of 
aliasing is to produce the Moire effect when the image is not properly sampled. 

It is also worth pointing out that the Shannon theorem states that the signal should 
not only be properly sampled, but also properly integrated over the sampling interval. 
In other words, one should have a 100% fill cycle (the fill cycle being the ratio of the 
integration time to the sampling time). If this is not the case, frequencies from the non- 
band limited signal above the Nyquist frequency will not be suppressed. It is then very 
critical to make sure that the fill cycle is as high as possible, thereby helping to reduce 
the aliasing effects of the high-frequency power contained in the non-band limited 
signal. The aliasing effects will severely affect the derivation of how the power drops 
at high frequency, implying that extreme care should be taken when recovering the 
power law of the background at high frequencies. From this point of vie w, the VIRGO 



instr ument aboard SOHO has, by d esign, a fill cycle ver y close to 100% (jFrohlich et al.l 
Il997h . while that of GOLF is 80% (|Gabriel et al.lll995l ). 

Last, but not least, since the signal is digitised, it is possible to use digital filters 
that are routinely us ed in the electro nic world for many digital applications (telephone, 
images and so forth; lAntonioulll979h . 

4-2.2 Spatial filters: Sun-as-a-star observations 

When the Sun is observed as a star, there is a cancellation effect due to the inte- 
gration over the disk. The net effect is that not all modes can be detected. The 
calculation of the attenuation a s a fu nction of the degree of the mod es is given by 
C hri stensen-Dalsgaard fc Goughl <|l982h for solar radial velocities and bv lToutain et al] 



(| 1999( 1 'or intensity; mainly modes with I < 4 can be detected. To the first order, when 
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the solar inclination angle is 90 degrees with respect to the line of sight, only the 
odd / + m modes are detectable. Neve rtheless, there is a smal l effect with the solar 
inclination angle discussed recently by iGizon fc Solanki d200ot). render ing all 21 + 1 



singlet modes of a multiplet detectable (See also lBallot et al, I l2006l , l2008h . Indeed 



sev- 



eral attempts had already been made to determine from the relative amplitudes of the 
observed p-mode singlets in the Sun t he orientation of the principal axis of rotation 
( Gough et al1ll995l ; | Gizon et al.lll998f) For g modes, th e sensitivity to the degree of 



the modes was derived by iBerthomieu fc Provost! |l990f ) for both intensity and solar 



radial velocities. They showed that for g modes, sensitivities of different degrees de- 
pends strongly on frequency. The large difference in sensitivities between solar radial 
velocities and intensities can in principle be exploited to help identify detected g modes. 

4-2.3 Spatial filters: spatially resolved solar images 

The observation of resolved images of the Sun makes possible the decomposition of 
the spatial structure of the modes onto functions that are related to the modes. The 
spatial filters may be classified in terms of different types of mask: 

— spherical harmonics masks 

— g-mode specific masks 

— optimal masks 

— time-distance masks 

For p modes, the obvious choice for the masks is to use the spherical harmonics that 
match the spatial dependence of the modes when the star is spherically symmetric 



.spatial 

( Unno et al] 1 19891 ). When the star rotates st eadily or has a magnetic field, the spa- 



tial dependence becomes more complicated ( Gough fc Kosovichevl [l993l ; IGizon et al.l 



19981 ) partly because the star is no longer spherically symmetrical, and partly because 
the modes must be referred to a frame of reference, if one exists, in which the state 
of the star is steady. If no such state exists, a well defined frequency of oscillations 
does not exist either. One has to bear in mind that such effects could be important 
for the visibilities of g modes if there is a la rge magnetic field in the central core 
( Goode fc Thompsonlll992l ; iRashba et al.ll2006n . 



Spherical-harmonic masks correspond to displacement or temperature perturba- 
tio ns (velocity or intensity observations). These masks were applied first in the 1980s 
by iBrownl \ 19851 ) (using an algorithm based on the Fast Fourier transform). The de- 



gree s ensiti vities of these masks h ave been derived bv lSchoul 1 19921 ) , lAppourch aux et al 



<|2000r i and ICorbard et all (|2008h for imaging instruments such as SOI/MDI, GONG, 
LOI, or SODIS]v(3 on the PICARD mission. 

Unfortunately, although the spherical-harmonic masks are rather well adapted for 
most p modes, they are not well adapted for the g modes. For velocity observations, 
the significant horizontal g-mode displacement must be taken into account (the per- 
turbations induced by p modes are predominantly radial), the spatial dependence of 
which differs from that of the radial displacement. For intensity observations, the os- 
cillat ions perturb the figure of the Sun modulating the light emitted from the sur- 
face (|Berthomieu fc Provost]|l990l ; iToutain et al.lll999h . For either t ype of observation , 



the additional contribution is a non-spherical harmonic function ( Unno et al.lll989l : 



lAppourchaux fc Andersenlll990l ). Specific g-mode masks can be devised (for intensity 
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observations see lAppourchaux fc Andersenlll990h . In some cases inclusion of the most 
appropriate spatial dependence of the g-mode perturbation still remains to be done. A 
large uncertainty lies in the way nonadiabatic effects are accounted for. 

Masks that are optimal for som e specific signals have been developed. For instance, 
lAppourchaux fc AndersenI ( 1 99Cf ) developed p-mode intensity masks that minimise 
both the leakage from other degrees and also the leakage from other m (see also 



Christen sen-Dal sgaardl 1 19841: Kosovichevj 



1986a 



1984T). Appourchaux et al l|l998h and by iToutain fc Kosovichevl 



Similar masks have been derived 



bvlGoueh fc Latour | 

(|2000h . Recently. I Wachter et all (|2002T l produced g-mode masks optimising (i.e. min- 
imising) the noise contribution from the supergranulation noise across the solar disc 
(which has a strong horizontal component). 

Other kinds of masks derived from the properties of the observations themselves 
have been used by IVecchio et al.l (120051 ). These masks are based upon a Proper Or- 
thogonal Decomposition (POD) of the velocity field, which is nothing less than the 
computation of eigenvalues and eigenvectors of that vector field. This POD analysis 
showed different power spectra at disc centre compared to the solar limb. These dif- 
ferences could be related to the dependence of the perturbation with distance to the 
centre of the solar disc. 

The aforementioned masks are defined for obser ving perturbations at the surface. 
With the advent of time-distance helioseismology ( Duvall et alj fl993h . there is the 
possibility of constructing special masks that will be sensitive to perturbations located 
deeper in the Sun, where the motions induced by gravity modes are much larger than 
at the surface. The idea is to detect motions induced by the g modes on the solar p 
modes. 



4.3 Spectrum estimators 
4-3.1 For data with no gaps 

The estimation of the spectrum of the time series is of prime importance when one 
wants to detect eigenmodes. The lifetimes of the modes should be taken into account 
when using these estimators. The following estimators are at our disposal: 

— Fourier spectrum (power spectra) 

— Lomb-Scargle periodogram and sine wave fitting 

— Average, smoothed and multitapered spectra 

— Cross spectrum 

— Random Lag Singular (Cross) Spectrum analysis 

— Frequency matching (oversampling and bin shifting) 

— Time-frequency spectrum 

— Varying time base 

Fourier spectrum estimation is widely used in helioseismolo gy. Its properties are 
well known and quite of ten well understood (|Bracewelll bOQOh , as are its statistics 
( Davenport fc RootHl958T ). In short, the Discrete Fourier spectrum estimation is widely 
used for time series being sampled at a regular sampling cadence (At) during an ob- 
serving time (T) that is an integer number TV of the sampling cadence. The Nyquist 
frequency of the Discrete Fourier spectrum is then i^N yqu i s t = 1 /^it , and the so-called 
frequency bin is the iVth part of that or Av = 1/N/At. The increase in frequency 
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resolution makes the detection of sine waves extremely effective because the power 
spectrum level of any stationnary noise decreases like 1/T. If a sine wave would change 
many times its phase during the observation time T, it would not be advised to com- 
pute the Fourier transform of the whole observation. In that latter case, the optimal 
computation would be to perform the Fourier transform for the period of constant 
phase, if these are known, and then to add the power spectra. An other incentive for 
not computing the whole Fourier transform would be be that the frequency of the sine 
wave changes gradually with time (See later in this Section). 

Other tools have been developed for time series that are not evenly sampled. Then, 
the Fast Fourier Transform should be replaced by an adapt ation of the D iscrete Fourier 
Transform, the so-called Lomb-Scargle (LS) periodogram |Scarglelll982r >. The LS peri- 
odogram is widely used in astrophysics for stars having modes with very long lifetimes. 
The LS periodogra m is strictly equivalent to sine-wave fitting, as shown in Appendix 
C of |S 3 (|l982T ). The main difference between the two methods is more in the ap- 



proach than in the result: with the LS periodogram the frequency of the sine wave is 
found from the spectrum itself, while for the sine wave fit the frequency is obtained 
from the fit provided that a proper initial guess is given, coming from the LS peri- 
odogram for instance. For spe ed, the LS periodogram is usually computed using the 



implementation prescribed by IPress fc Rvbickil ( 1989T ) which is an approximation of 
the LS periodogram based upon extirvolation 1 ™ on a regular mesh and the use of the 
FFT. The prescription is then very close to interpolating onto a regular mesh. It is 
worth noting that most scientists using the LS periodogram for unevenly sampled data 
are in fact computing the FFT of the original data resample d onto a regular mesh, but 



with only a proper normalization as given bv lScarglel l| 19821 ) 



Fourier spectrum estimation is well adapted for periodic signals (pure sine waves 
or stochastic waves) but not necessarily well suited for estimating the spectral density 
of pink or red noise. For that purpose, one can: 

— average the power spectrum over an ensemble of q' sub-series, 

— smooth the power spectra over q' frequency bins or, 

— use multitapered spectra using the full time series for deriving a similar average. 

Fourier spectrum estimation is now bein g replaced by m ultitapered spectra that are 
widely used in geophysics (for a review see lThomsonlll982T ). Multitapered spectra are 
generated by applying a set of slepian tapers to a single time series, and an estimate of 
the mean power spectrum is derived from an average of these spectra. The multitapered 
spectra are statistically independent from one another, and the statistics of the mean 
spectrum follows a y 2 dist ribution with 2q degrees of freedom (where q is the number 
of tapers, I Thomson! \l98$ ). While the statistics of the average power spectrum (or 
smoothed power spectrum) also follow a \ 2 distribution with 2q degrees of freedom, 
the resolution of the average spectrum is q' times lower than that of the multitapered 
spectrum. In helioseismology the use of t hese slepian tapers h as been replaced by more 
practical (but less accurate) sine tapers ( Komm et al.lll999l ) . Unfortunately, for long- 



lived modes, tapers tend to broaden the peaks, as shown by iThomsonl ( 19821 ). Tapers 



as such provide more benefit for p modes having a lifetime shorter than the observation 
time. 

The use of cross spectrum esti mation can also be us eful for detecting signals being 
present in two different time series (|Sturrock et edj|2005h ■ This can be done either from 



17 Reverse interpolation or extirpolation replaces a function value at any arbitrary point by 
several function values on a regular mesh 
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the same observable or from dif ferent observables (e .g. intensity and radial velocity). 
This technique has been used bv lGarci'a et al.1 (jl999h for improving the signal-to-noise 



ratio in the p-mod e regio n, but it is not useful for g-mode detection as shown by 



lAppourchaux et al, I (|2007T ) 



Random Lag Singular (Cross) Spectrum analys is is an elaborate techniq ue based on 
Singular Value Decomposition (RLSSA, RLSCSA Ivaradi et al.lfl999l . l200oh . Although 



the technique claimed suc cessful detection of low-frequency p modes (jBertello et al] 
bOQObl ; [Garcia et al.lliooTbh . although a proper assessment of its statistical properties 
is still lackin g: i.e. it cannot be excluded that the technique produces large peaks solely 
due to noise ijCouvidatlEooi ). The spectrum produced by the RLSSA technique can be 



compared to a power spectrum raised to an unknown power (greater than one). As such, 
small peaks with a signal-to-noise ratio slightly larger than unity are amplified, while 
those with a lower signal-to-noise ratio are damped. The fact that this exponentiation 
is unknown (dependent upon the data) makes the RLSSA and RLSCSA techniques in 

our view unsuited for robust signal analysis. 

Frequency matching has been developed bv lGabriel et~aT] |2002l ) and bv lChapl in et al.l 



( 20021 ) using zero padding and frequency-bin shifting, respectively. When one wants to 



detect signals from oscillators having a lifetime longer than the observing time, there is 
a significant chance that the frequency bin will not match the frequency of the oscillator. 
As a result, power from the signal may be split between f requency bins and th e main 
peak could be reduced in power by up to 60 %, as shown bv lGabriel et al.l ( 20021 ) . In or- 



der to alleviate this problem one can either oversample the data by using zero padding 
Gabriel et aHl2002h or try to tune for the frequency by using frequency-bin shifting 



Chaplin et alil2002l ). This latter technique involves creating many dif ferent time 



ries o f similar but slightly different lengths, but without zero padding (IChaplin et al.l 
I2002T ). In either case, the statistics of the observation are indeed affected. For oversam- 
pling, the analytical calculation has been replaced by Monte-Carlo simulation showing 
that there are typically three indepen dent frequency bins w hen the Fourier spectrum is 
oversampled by a factor five or more (jGabriel et al.ll2002h . The frequency-bin shifting 



method also produces three independent power spectra from the many spectra gen- 
erated (|ChapIin et al.ll2002h . Both approaches have the added benefit that they give 



access to different "realizations" of the noise background, as if we had three actual, 

independent realizations of the noise available. 

Time-frequency spectrum o r wavel et analysis has been used b v lGabriel et al.l ( 19981 ) , 



bv lFinsterle fc FrohlicrJ (|200ll ) and by I Jimenez fc Garcfal l|2009l ) to look for signatures 



that are affected by the lifetimes of candidate g modes. The statistical properties of 
these wavelets can be derived from those of the Fourier spectrum. The potential influ- 
ence of solar activity on the size of the resonant cavity, which w ould affect the frequency 
of long-lived modes, has also been studied by iGabriell (|2006l ). If the frequency of the 



mode changes slowly with time, the freq uency of the m ode will be spread over sev- 
eral frequency bins. The idea followed by iGabriel ( 20061 ) is to change, or distort, the 



sampling time to follow the variations in the mode frequency induced by solar activity. 
4-3.2 For data with gaps 

All the estimators described above have different properties when gaps appear in the 
data. Hereafter, we shall discuss mainly the following: 



— Fourier spectrum (power spectra) 
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— Lomb-Scargle periodogram and sine wave fitting 

— Average, smoothed and multitapered spectra 

— Frequency matching (oversampling and bin shifting) 

The impact of the gaps on the Fourier spectrum has been described bv lGabriel (|l994l ^ 
providing as a result that the correlation matrix of the frequency bins is not longer 
diagonal. The gaps introduce correlation between the frequency bins that require to 
be taken int o account when one w ants, for example, to fit the power spectrum (See for 
applications Istahn fc Gizonl [20081 ). 

The LS periodogram does not provide a better solution to coping with the presence 
of gaps. The reason is that although the Fourier transform explicitly includes gaps 
as zeros, adding zeros is also implicitly performed with the LS periodogram. As a 
consequence, correlations between frequency bins also exist with the LS periodogram, 
but these are frequently ignored. 

The impact of the gaps for the average will have the same origin as for the Fourier 
spectra. As for the smoothed spectra, there is an intrinsic correlation between the q' 
bins; the number of truly independent frequency bins is then divided by q' . In that 
latter case, the gaps do not have a large influence on the smoothed spectrum unless 
the fraction of gaps is about 1/q . For multitapered spectra, the influence of gaps can 
be ta ken into account for obtaining optimised tapers that match the structure of the 
gaps (jFodor fc Starklll99Sl ). In that case, the correlation between the frequency bins, 



although reduced, will be not negligible, especially if several tapers are used for the 
estimation of the mean power spectrum. 



4.4 Statistical testing 

Statistical testing is essential when one wants to decide: have we found g modes or 
not? This is related to decision theory, which can be summarised as how do we choose 
between one hypothesis versus another in the presence of uncertainties? In this area, 
there are two schools of thought: the frequentist school and the Bayesian school. 

The difference between a Bayesian and a frequentist relates to his or her views of 
subjective versus objective probabilities. A frequentist thinks that the laws of physics 
are deterministic, while a Bayesian ascribes a belief that the laws of physi cs are true 
or ope rational. The subjective approach to probability was first coined by |Pe Finettl 
(| 19371 ). For the rest of us, the difference in views can be summarised by this quote from 
the Wikipedia encyclopaedia: Whereas a frequentist and a Bayesian might both assign 
a probability ^ to the event of getting a head when a coin is tossed, only a Bayesian 
might assign a probability jqqq to personal belief in the proposition that there was life 
on Mars a billion years ago, without intending to assert anything about any relative 
frequency. In short, frequentists assign probability to measurable events that can be 
measured an infinite number of times, while Bayesians assign probability to events that 
cannot be measured, like the outcom e of sport-related bets, for instance, or the survival 
time of the human race |Gottll 19941 ). 



With this word of caution, one should never forget that we wish to use this armada 
of statistics because we know that we are at the limit of detection for g modes. In what 
follows, we will try to give an overview of what we believe we know on a subject that 
is fast evolving; and what we write is certainly not gospel. 
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4-4-1 Frequentist approach 

For a frequentist, statistical testing is related to hypothesis testing. In short, we have 
two types of hypotheses: 

— Ho hypothesis or null hypothesis: what has been observed is pure noise? 

— Hi hypothesis or alternative hypothesis: what has been observed is a signal? 

For the Hq hypothesis, we assume a known statistics for the random variable X ob- 
served as x and assumed to be pure noise; and then set a false ala rm probability that 



defines the acceptance or rejection of the hypothesis (|Scarglelll98l ). The so-called de- 
tection significance (or p-value, terms not widely used in helioseismology) is the prob- 
ability of having a value as extreme as the one actually observed. There is an on-going 
confusion because statisticians call the significance level what astronomers call the false 
alarm probability; and they call the p-value what is set in astronomy as the detection 
significance (which is not the significance level). Here we shall use the current vocab- 
ulary understood in astronomy. For example, the false alarm probability p for the Hq 
hypothesis is defined as: 

p = P (T(X)>T(x c )), (29) 

where T is the statistical test, and Po is the probability of having T(X) > T(x^) when 
Ho is true; and x c is the cut-off threshold derived from the test T and the value p. 
For example, take the case of a random variable X distributed with x 2 , 2 degrees of 
freedom (d.o.f) statistics, having a mean of a. If we further assume that T(X) = X, 
we then have that: 

P= Po(X >x c ) =e~^. (30) 

If one observes a value x of the random variable X that is larger than x c , the Ho 
hypothesis is rejected. The value that is quoted in this case is the detection significance 
V, i.e., 

V = e~°. (31) 
The Ho hypothesis was used by lAppourchaux et al. to impose an upper limit 



on the g-mode amplitudes. The method was based on the knowledge of the statistical 
distribution of the power spectrum of full-disc integrated instruments, namely the x 2 
distribution with 2 d.o.f. 

For the Hi hypothesis, we assume given statistics both for the noise and for the 
signal that we wish to detect, and set a level that defines the acceptan ce or rejection 
of that hypothesis. The Hi hypothesis was used in iGabriel et al] |2002l ) to determine 
the p robability of detecting a sine wave given noise in the GOLF data ( Gabriel et alj 
I2002T ). In this latter case, the statistics follow a non-centred \ 2 distribution with 2 
d.o.f. Both hypo theses have also bee n used to calculate significance levels for detecting 
stellar p modes ( Appourchauxll2"o04l ) . 



Taking a decision based on the result given by a single test, for either hypothesis, 
could lead to errors in the decision process. For instance, the null hypothesis could 
be wrongly rejected when it is true (false positive or wrong detection), but could also 
be wrongly accepted while it is false (false negative or no detection in presence of a 
signal). The false positive results in a Type / error, while the false negative results in a 
Type II error (see Table 2) . The ideal case would be to set a test that would minimise 
the occurrence of both types of errors. 

In helioseismology, it has been cu stomary when applying th e Hq hypothesis to set 



the decision level arbitrarily at 10% ( Appourchaux et al.ll200oT ). From the frequentist 
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Table 2 Types of error obtained for different decisions, based upon the statistical test per- 
formed, and how the error relates to the status of the Hq hypothesis. 





Status of Hq 




True 


False 


Reject 


Type I 


Correct 


Decision 






Accept 


Correct 


Type II 



view point there is nothing wrong in setting a priori the decision level before the test 
is applied. There are three types of result we might obtain from the test: 

1. Ho always rejected 

2. Ho rejected or accepted at a level very close to 10% 

3. Ho always accepted 

Decision 1 will lead to the mention of a detection being statistically significant at a 
level provided by the detection significance [from Eq. (|31|l ]. The next step would be 
the application of a test for the Hi hypothesis, very likely resulting in the detection 
of signal. Decision 2 is the more difficult borderline case, forcing us to either accept or 
reject Hq. Here, we might ask: are things really that clear-cut? What are the chances 
that if we accept Hq is it actually wrong (Type II error), or truly right if rejected (Type 
I error)? Decision 3 seems straightforward, i.e., noise dominates, but might one then 
be tempted to lower, a posteriori, the decision level? 

These potential actions result from the application of a frequentist test trying to 
answer the following question: what is the likelihood of the observed data set x, given 
that Ho is true or p(3:|Ho)? The detection significance mentioned when the test rejects 
the Ho hypothesis is nothing but p(a;|Ho), when actually what we want to know is the 
likelihood that Ho is true given the data, i.e., p(Ho|x) T>). The frequentist view 
does provide a useful answer when one can repeat the observations ad infinitum. But 
when we have only one observations, another approach we may apply is the Bayesian 
approach, which in principle gives access directly to p(Ho|x). 



4-4-2 Bayesian approach 



On the posterior probability We should never forget the two sides of the coin: if 
probability (likelihood) can justify alone the rejection or acceptance of an hypothesis, 
this probability is not the significance that the hypothesis is rejected or accepted. 
The decision levels discussed above are related directly to a well-known controversy 
in the medical field, concerning improper use of Fi sher's p-values as measures of the 
probability of effectiveness of a medicine or drug ( Sellke et alj l200lh . The detection 



significance (or p-value) is improperly used as the significance of the evidence against 
the null hypothesis. It is far from trivial at first sight to understand what is wrong with 
the detection significance. Let us recall the example of a noise with 2 d.o.f statistics. 
In that case the detection significance is given as: 

V = e~° j£ P Q {X > St). (32) 



The latter statement (^) is fundamental. The observation is performed only once pro- 
viding a value of Sc. It is not correct to assume that this observation if it were repeated 
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would provide the same level x. The mistake is to ascribe a significance to a measure- 
ment performed only once, i.e., not repeated, and spanning just a very small volume 
of the space of the parameters (e.g. X £ [x, x + 5x]). If one makes a measurement x of 
the random variable X which is above x c , the significance of that measurement is not 
g-a/cr j n ^ e f ramewor k f Bayesian statistics, we are not interested in the detection 
significance but in the posterior probability of the hypothesis, in other words as already 
stated above p(Ho|ai) ^ T>. A similar description of this misunderstanding has been 
presented bv lsturrock fc Scarglel |2009h . 



In order to derive the posterior probability p(Hg|a:), let us first recall the Bayes 
theorem. The posterior probability of a hypothesis H, given the data D and all other 
prior information I is stated as: 

(33) 

where P(H|I) is the prior probability of H given I, or otherwise known as the prior; 
P(D|I) is the probability of the data given I, which is usually taken as a normalis- 
ing constant; P ( H|I) i s the direct probability of obtaining the data given H and I. 
IBerger fc SeUkel jl987t ) obtained, using the Bayes theorem, p(Hg|i) with respect to 



p(i|Ho) and p(x|Hi), where Hi is the alternative hypothesis. 



/tt|-\ p(H )p(a|H ) 

P[no}X> p(Ho)p(i|H )+ P (H 1 )p(£|H 1 )- (M) 



We set po — p(Hq), and since we have p(Hi) = 1 — po, they finally obtained 

(1-po 



(35) 

V Po J 

with C being the likelihood ratio defined as: 

£ = Pj 0§ 1 y (36) 
p(x|Ho) 

Here, p(Ho|i) is the so-called posterior probability of Ho given the observed data St. 
Naturally there is no way to privilege Ho over Hi, or vice versa, otherwise our own 
prejudice would m ost likely be confirmed by the test, i.e. po = 0.5. Subsequently, 
IBerger et all l| 19971 ) recommended to report the following when performing hypothesis 



testing: 

if £ < 1, reject Ho and report p(Hrj|i) = T~~£ i (37) 

if C > 1, accept Ho and report p(Hi|i) = - ■ (38) 

The advantage of such a presentation is that even for a borderline case, say when the 
ratios above are close to unity, it is clear that there is only a 50 % chance that the Ho 
hypothesis is wrongly accepted, or wrongly rejected. This presentation is more honest 
and better enca psulates hum a n jud gement and prejudice. 

The work of lBerger et"al] (|l997l ) can be applied to the problem of detecting modes 



with very long lifetimes, i.e., those modes whose peaks are predominantly restricted to 
one frequency bin in the frequency spectrum. We can write the likelihood for the Ho 
hypothesis as: 

P(S|Ho) = \^' B - (39) 
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where B is the mean noise level. Next, we assume that the mode is stochastically 
excited with a known amplitude A. The likelihood for Hi is then: 

p(2|Hi) = bTa^' (B+A) - (40) 

Since the detection significance is p = e~ x l B , we use Eq. 1)350 to give: 

-1 



p(Ho|x-)=(l + ^p- A/(B ^M 



(41) 



Figure [16] show the results for two different detection significances. When the detection 
significance is 10 %, the likelihood ratio can be greater than unity for large values of 
the mode amplitude, leading to the acceptance of the null hypothesis. This is rather 
paradoxical, i.e., that large mode amplitude can lead to the rejection of the alternative 
hypothesis. To resolve the paradox we note that the posterior probability of Ho is in any 
case never lower than 40%, or the posterior probability of Hi is never higher than 60%. 
This implies that both hypotheses are equally likely when the detection significance is 
as low as 10 %. In other words, when we set, a priori, a large mode amplitude and get a 
low detection significance, the alternative hypothesis is as likely as the null hypothesis. 

The main conclusion to be drawn from this calculation is that the detection signif- 
icance should be set much lower than 10 % in order to avoid misinterpretation of the 
result. For example, with a detection significance of 1 %, the post erior probabi l ity for 
Ho can fall to 10% when the signal-to-noise ratio is above unity. fSellke et all ((2001) 
showed that the posterior probability can never be lower than the lower bound 



^H 1 -^; • (42) 

The reader may verify for themselves that this lower bound is effectively reached for 
Eq. (|41|l . In the case, when the amplitude of the mode A is not known, one needs to 
set, a priori, value for the likely range of amplitude. In the case of a uniform prior, the 
posterior probability p(Ho|:r) then does reach a minimum that is higher than the lower 
bound of Eq. (|42jl . 

In summary, the significance level should not be used for justi fying a detection (o r 
a non-detection). Instead we recommend to use the prescription of iBereer et all l|l997l ) , 
as given by Eqs. 1)370 and (|38)l and to specify the alternative hypothesis Hi. 



On the choice of the prior probability One important question when applying 
Bayesian statistics is what value should the prior probability of the hypothesis Ho, i.e., 
Po, take? We define the prior probability as the probability that the Hq hypothesis 
is correct. The probability that the alternative Hi hypothesis is correct can then be 
defined as p(Hi) = 1 — p$. It is common to set po = 0.5 so as to avoid prejudicing one 
hypothesis over the other. Given current model prejudices, we kn ow that low-frequenc y 
p modes have narrow widths typically between 10 and 30 nHz, (jchaplin et alJl2005l ). 



which cover a few frequency bins only; and we expect the same to be true for g modes. 
In which case, would we expect the probability that Hi and Ho are true to be the same 
at all frequencies? Since Bayesian statistics uses a priori knowledge, it is possible to use 
our knowledge of the properties of modes to tell us which hypothesis is more likely to 
be true at a given frequency. We can also use this knowledge to then guide the regions 
we search in frequency (See section 4.5.3). 
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Mode amplitude 



Mode amplitude 



Fig. 16 On the left-hand side, likelihood ratio £ as a function of the mode amplitude for 
detection significances of 10% (solid line), and of 1 % (dashed line); the noise is set to unity. 
On the right-hand side, the posterior probability of Ho as a function of mode amplitude for 
detection significances of 10% (solid line), and of 1 % (dashed line) (from Eq. I41H 1. 



4.5 Patterns 

4-5.1 Rotational splitting 

As shown in Section 2 [see Eq. (I19|l ]. the mode degeneracy is lifted by the solar rotation; 
each (n, I) mode is then split into 21 + 1 components. The collapsogramme technique pio- 
neered bv lAppourchaux et alJ 1 2000l ) makes use of the pattern created by the rotati onal 



splitting to detect the modes. It has been extensively used bv lSalabert et al.1 (|2009h for 



detecting modes below 1000 /iHz, using resolved-Sun data collected the GONG instru- 
ment. It can also be used for full-disc instruments producing a sing le power spectrum: 
in th is case it is called an overlapoqramme (For an application see ISalabert fc Garcial 
120081 ). IChaplin et al.1 (|2002T l has devised a statistical technique based on the detection 
of an ordered multiplet (due to rotational splitting) in the power spectrum of full-disc 
integrated data. This lowers the detection level, relative to the level for just a single 
peak, depending on the number of peaks that are found in the multiplet. The derived 
limit (under the Ho hypothesis) can be translated into the equivalent a level in a power 
spectrum for a one-year observing time (where a is the mean power spectral density in 
the neighbourhood of the sought-for modes): 5.9<r, 4.5a, 3.8a for a doublet, tri p let and 
a quadruplet, respectively. A similar approach has been used bv lGarcla et al. (2001a) 



for the GOLF data. They derived levels for detecting modes of various degrees using 
Monte-Carlo simulations, which can be translated for a one-year time series to 5.4a for 
an I = 1 doublet, 5.9a for an I = 2 doublet and 4.3a for an I — 2 triplet. 

4-5.2 Asymptotic behaviour 

The asymptotic behaviour of g-mode periods may be used to aid detection [see Eq. 
of Section 2]. Such a nalyses were pioneered by Delache in 1983, leading to a claimed 



detection of g modes ( Delache fc Scherrerlll983T ). This approach is only of relevance to 



high-order g modes (i.e., very low-frequency modes below ~ 100 /iHz) for which the 
asymptotic behaviour applies. Unfortunately, the solar noise increases towards lower 
frequencies, and the mode spacing (in frequency) decreases dramatically. The situation 
is further complicated by the effects of rotational splitting (|Frohlich fc Andersen|[l995l . 
see also Section 5.2.1). 
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The E xact Fractions Tech nique (EFT) was pioneered by Ivan der Ra av (198^) and 



Palle et "all ( 19981 ) . The technique is rather similar to the echelle diagram 



Gred jl98ll ). The main difference is that instead of cutting the spectrum 



applied by 
devised by 

into several windows of width Avq, the frequency axis is first changed to period, and 
the spectrum is then cut into several windows of width P p/\/l(l + 1). The full EFT 
procedure is described in more detail in lPalle et al.l ( 19981 ). 



The EFT uses the property that the modes , if they exis t , are regularly spaced 
in period. The same property has been used by iGarci'a et al. (2007). They compute 



the periodogram of the frequency power spectrum expressed in period, i.e., the pe- 
riodogram of the periodogram. The regular comb-like pattern expected in the first 
periodogram then manifests as a nother comb-like pattern i n the second periodogram. 
Unlike the technique devised by iDelache fc Scherrerl ( 198&t ) or the EFT, the statisti- 



cal understanding is far from trivial, requiring Monte-Carlo simulations to derive the 
likelihood of structures appearing above a given threshold. 

4-5.3 Guided search 

An artificial way of reducing the detection limit is to reduce the size of the frequency 
windows over which we wish to search f or modes, e.g., by looking in windows centred 
around theoretical g-mode frequencies. iDenison fc Waldenl (| 19991 ) provided a simple 
formula to derive the number of peaks due to noise that one can find in a power 
spectrum, given a list of frequencies and a window containing these frequencies. Under 
the Hq, hypothesis, it is written as follows: 

N = N l (l-(l-p det ) N "). (43) 

Here, N\ is the number of frequencies guiding the search, p^ ct is the probability level 
needed for identifying a peak and N w is the window size in units of frequency bins. 
When NwPdet is much smaller than unity, we can rewrite Eq (|43p as: 



N « A/iJVwPdet- (44) 

This simple formula is quite useful, in that it allows us to realise that the number of 
identified peaks will increase with the size of the window and the number of frequencies 
guiding the search. This is the drawback of such a method: spurious peaks will be 
detected that are likely to be wrongly identified as g modes. Here we should also 
remind the reader that theoretical p-mode frequencies showed systematic errors of the 
order of several /xHz, which came from our inability to model properly the surface of 
the Sun ( Christensen- Dalsgaardlfl990l) . Therefore care should be exercised when using 



theoretical frequencies as a guideline for searching for g modes. 

The sensitivity of g-mode frequency predictions to solar models is of the order of 
1% (see Section 2.2.3). This error estimate can be used to produce an inverted top- 
hat prior probability distribution that can be used to guide Bayesian searches for g 
modes. The cent res of each hat were d etermined by the mode frequencies predicted by 
the Ml model of IProvost et al. I (|2000l) and the width of each hat was 2% of the mode 



frequency. Frequencies where modes are expected to be observed have been given a prior 
probability that the Ho hypothesis is true (po m equation 36) of 0.5 and elsewhere the 
prior probability is unity. Figure [17] shows such a distribution, where all components 
that will be visible in Sun-as-a-star data (when I + m is even) have been plotted for 
I — 1 and I = 2 g modes. We have assumed that the different m components are 
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Fig. 17 An Echelle diagram showing the variation of the of the prior probability for g modes. 
The prior probability, prji wa s either 1.0, at frequencies where modes are not expected to be 
observed, or 0.5 at frequencies where modes are expected. 



separated by 0.4 piHz. However, as we do not know the rotation rate of the solar core 
this assumption may not be valid (see Section 2.2.4). Notice that despite the large 
number of g modes that are present at low frequencies the prior is not uniform with 
frequency and so this approach should limit the number of false detections. 



4.6 Data combination 

The signal from solar oscillations will be common to contemporaneous data observed 
by different instruments. It is, therefore, pertinent to search spectra constructed from 
contemporaneous data for statistically significant prominent concentrations of power, 
which lie significantly above the local noise background and that are coincident in 
frequency in the different datasets. The use of different data sets related to different 
observables and/or wavelengths could very well be the solution to the g-mode detection. 
There is no doubt that the combination of more than two signals could considerably 
lower our detection limit. Observables such as radial velocity, intensity fluctuations, 
limb displacement and/or brightening are polluted by different sources of noise, such 
as supergranulation and active regions, that produce different signatures. We can list 
the possible combinations as follows: 

— use of one instrument 

— use of more than one instrument 

In terms of the analysis techniques we might apply to exploit such data, autocorrela- 
tion fa lls into the first category. The collapsogramme technique used bv lAppourchaux et alj 
(20001) also belongs to the first category. Another technique developed bv lGarcia et al. 
( 1999 ) used a longer sampling time (e.g., 80 sec instead of 40 sec) to create two over- 
lapping, independent timeseries of the GOLF data. The first time series used all the 
even-sample data points, while the other used the odd-sample data points. The combi- 
nation of the two resulting Fourier spectra can be used to improve the signal-to-noise 
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ratio of modes with a lifetime shorter than the observation time ( Appourchaux et al.l 
120071 ). 

Th e Multivariate Spectral Regression Analysi s (MSRA) belongs to the second cat- 
egory 1 Koopmanslll974 : lAppourchaux et al]l200oh . It assumes that the modes are pre- 
dominantly coherent over the observing time. The basic assumption is that low-degree 
p modes and g modes will have lifetimes that are significantly longer than the observ- 
ing time. The M SRA has been applied to da ta from the different SPM channels of the 
VIRGO data bv lFinsterle fc Frohlichl (|200lh . 

A derivation of the joint prob ability, p, of the o ccurr ence of coincident prominent 
features i n two spectra is given in lBroomhall et aL Various statistical tests are 

derived in IChapl in et al.l (|2002h that determine the probability of observing prominent 
structures in a single spectrum. These tests can be easily adapted so they can be used 
to compare two spectra if the probability of observing a si ngle prominent spike a t 
the same frequency in each of the spectra can be calculated ( Broomhall et alj|2007l ). 
The technique works best when the signal-to-noise ratio is similar in the different 
instruments. These tests were originally designed to detect p modes, however, the 
underlying theory is still applicable when searching for g modes. 



5 Applications of detection techniques 

As outlined in the previous section, the building of a detection is made of an assemblage 
of various elementary bricks. It is essential that a detection technique be segmented in 
different categories. The key category is statistical testing because it exemplifies that 
the potential biases of a detection technique can indeed be understood. 

Hereafter, we will try to review the application of detection techniques, and their 
associated results. 



5.1 Statistical thresholding 
5.1.1 Radial velocity 



The first reported g-mode detection was made by iDelache fc Scherrerl 1 19831 ). They 



analysed differences in the solar radial velocity recorded at disc centre, and further 
out toward the limb . Observations were made at the Stanford Sol ar Observatory 
dScherrer et al.lll979l ). and at the Crimean Astrophysical Observatory (|Severnvi et al.l 



1976). The length of the time series was 10-") days, but with a duty cycle of only 10%. 
After computing the power spectrum, they constructed the cumulativ e power distri- 
bution in the frequency range [45-105] /xHz. IDelache fc Scherrerl l| 19831 ) calculated the 



cumulative distribution of the power spectrum in this range, and ascribed the strongest 
peaks as being due to g modes. 

The cut-off used to tag prominent peaks was 2.5 a which, given the assumed x 2 
with 2 d.o.f statistics, means that in the window considered about 5 ± 2.5 "detections" 
will have been due solely to noise (there are about 55 independent frequency bins in 
the window given the 10% duty cycle). Fourteen possible modes were identified, at 
least half of which are likely due to noise. In addition, they did not take into account 
that the power spectrum has a 1/u dependence implying that the mean varies with 
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frequency; this explains the departure observed from the expected statistical distribu- 
tion. The minimum amplitude of the detected peaks was 25 cms" 1 . Since this paper 
appeared, no independent confirmation of the results has appeared in the literature. 
Indeed, as we shall see we below, quoted upper limits on g-mode amplitudes, given 
by analyses of mod e rn da tasets, lie some two orders of magnitude below the claimed 
IDelache fe Scherreri (jl983h amplitudes. 



A similar approach was used by lAppourchaux et al.l ( 2000l ) for setting an upper 
limit to the amplitude of the g mode s, which lay significantly below the amplitudes 



g mode s 

claimed by IDelache fc Scherre r (1983). The canonical figure for Hg was given for a 



threshold of 10 % for a 100/^Hz window, providing a minimum posterior probability for 
Ho of 38% (See Section 4.2.2). As shown in Fig. 17, a threshold of 1 % would provide 
a lower posterior probability for Hq of typically 15 % for mode amplitude at a signal- 
to-noise ratio less than 10 (but greater th an unity) . The 1 % t h resho ld would increase 
by a factor two the upp er limit quoted by lAppourchaux et al] ( 2000l ) . 



Gabriel et all (|2002l ) also used a statistical thresholding approach. They used over 



sampling in an attempt to increase the detectability of long-lived signals. They reported 
the detection of three g-mode candidates below 290 /iHz with a significance level of 
4%, or a posterior probability for Ho of 25%. They could not tag the candidates as 
being g modes, and then set an upper limit to the amplitude of g modes, based upon 
the Hi hypothesis, of 7 mm s ~ . 

iTurck-Chieze et al. d2004Tl applie d thresholding for various kinds of features: for 
singlets (as Appourchaux et al.ll200ol ). for doublets and for triplets. The latter appli- 



cation of the thresholding is more akin to searching for pattern of g-mode peaks, which 
is described below. Several candidate detections were flagged at a threshold of 10 % 
(implying a lower bound on the posterior probability of 38%). But none appeared in 
all their test observation periods. 

5.1.2 Solar diameter 

Observations of the solar diameter have also been used in the search for g modes, e.g., in 
observations by the SCLERA instrument. The lowest amplitude detected corresponded 
to a radius change of 1 mas for p modes, and to 0.1 mas for g modes; assuming that 
these radius changes were real physical changes, they would have been e quivalent to ra - 
dial velocities of 50 km s _1 and 250 m s _1 , respectively. As outlined bv lBrownl(|l979l ). 



the observed oscillation signals were more likely due to bright ness effects than to phys- 
i cal ra dius changes. This latt er fact was also confirmed by lAppourchaux fc Toutainl 
\ 19981 ) and lToner et~ai1 \ 1999^ 1 using space-borne instrumentation. The amplification is 



due to a geometrical effect, which makes optical thickness ch anges dominate tempera- 
ture c hanges. The rms amplitudes of the modes detected by I Appourchaux fc Toutainl 
(Il998l ) were t ypica lly of the order of 8 ppm, corresponding to a radius change of 1 mas. 



Toner et al. I (|l999h showed that the oscillation signal increases as the distance to the 



limb decreases. This property will be used b y the SODISM instrument on the PICA RD 
mission to try to detect modes at the limb ( Dame et al.ll 19991 : ICorbard et al.ll2008l ). 



5.1.3 Solar wind data 



Thomson et al. (1995) used magnetic field and particle data collected by t he Ulysses 



mission. They used different spectral estimators (based on multi tapers; see I Thomson! 
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( 19821 )) and different statistical tests (based on the F-test), with a threshold of 5 



They claimed the detection of s everal g-mode frequencie s. How ever, these claimed 
detect ions were not confirmed bv lDenison fc Walde 3 |l999l ) or bv lHooge fc Rilevl 



( 19981 ) . who both argued that the chosen threshold was too w eak and gave rise to false 



detections. By making use of the work of lBerger et alJ ( 19971 ). we can now understand 



in retrospect that this threshold will have lead to a posterior probability of Ho being 
true of at least 29 % [from Eq. (|42[) ] : the applied F-test threshold was clearly not low 
enough. 



5.2 Searching for patterns of g-mode peaks 
5.2.1 Asymptotic properties of periods 

The fir st attempt to use the asym ptotic properties of the g-mode frequencies dates 



theasymp 

back to lFrohlich fc Delachel (|l984al rbl). They used the asymptotic formula [see Eq. pip]. 



corrected for rotational splitting, to compute g-mode frequencies as a function of Pq 
and Hq (S7q being the mean rotation as defined by Eq. [19}. They then calculated the 
coherence between the calculated g-mode frequencies and the observed power spectrum, 
and drew maps of coherence as a function of (Pp, Op ) . No detections were made. The 
sa me procedure was appl i ed bylFrohlich fc Andersen! ( 1995I - ) also without success. 



Frohlich fc Andersenl (|l995h concluded that either the g modes were not present 



or that the solar noise was too high. This method has also recently been extended 
to include also g-mode frequencies excluded by the normal asymptotic equation. The 
models shown in Table 1 cover a sufficient range of Po rendering an interpolation 
scheme as function of Pq possible. This allows to include g modes with I — 1 and I = 2 
down to n = 1. The results are not yet conclusive, possibly due to the precision of the 
frequencies and/or the interpolation scheme. Another reason may have been that the 
range for the rotation in the core was too small. One of the major drawbacks of the 
technique is that the computation of g-mode frequencies needs to be very accurate. 

5.2.2 Rotational splitting 



Turck-Chieze et alJ (|2004l ) applied statistical thresholding to search for patterns of 



peaks induced by rotational splitting. They derived, using Monte-Carlo simulation, 
pattern rejection thresholds for doublets and triplets at the 10 % level for a 20-/JHz 
window (p osterior proba b ility greater than 38% ); similar levels were analytically 
derived by IChaplin et ail (|2002T l for p modes. Given the appeara nce, and disappear- 



ance, in different observation periods of the various detected peaks. iTurck-Chieze et all 
(2004]) concluded that some patterns were due to noise, thereby supporting the pos- 
terior probabi l ity of Ho being greater than 38 %. The analysis was further pursued by 
iMathur et ail l|2007l ) who confirmed detection at a 2% level for a 10/xHz window (corre- 



sponding to a 4 % lev el for a 20 /iHz window, or posterior probability greater than 26%) 
IChaplin et ail (j2002T l also searched for similar patterns, reporti ng possible det e ctions 



of p modes as low as 700 pHz. A similar technique was used bv lSalabert et all (|2009l 'l 
to detect low-order p modes, of higher degree, in resolved-Sun observations made by 



18 The .F-test is used for comparing two random variables having a \ 2 distribution with 2 
d.o.f 
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the GONG network. They detected modes below 1000 fiRz, with typical lifetimes of 
the order of 3 years. 



5.2.3 Periodogram of Fourier space 



Garcia et all (|2007h found a ver y interesting peak i n the periodogram of the peri- 



odogram of GOLF data (See also[G arci'a et aLll2008ah . They claimed that this peak is 
due to the superposition of g modes that lie in the asymptotic region, where the sepa- 
ration in period is approximately constant. The likelihood that the peak is not due to 
noise is 99.5±0.13%, with a posterior for Ho greater than 3.4 %. The peak is positioned 
close to the period predicted by standard solar models. The peak can be produced in 
several ways, all related to clusters of peaks being equally spaced in period. For in- 
stance, the power of the g modes could be spread onto several adjacent peaks like the 
presence of a deep magnetic field or changes in the propagation cavity (a displacement 
of the tacho cline as a con s equen ce of the activity cycle). Other possibilities were also 
discussed bv lGarcfa ei~al1 (|2007f ). They used artificial data to show that, if this peak is 



evidence of g modes, the modes could have widths in a frequency-power spectrum that 
are commensurate with damping times of several months, with a possible signal-to-noise 
ratio not greater than five in power. If this is the case, and given that the data searched 
here are 3071 d in length, the width of these modes would be resolved in the spectrum. 
For example, if the mode damping time was 122 d the power of the mo des in the spec- 
trum would be spread across eight frequency bins. lAppourchauxl ( 20041 ) showed that it 



is possible to detect short-lived modes by smoothing the power spectrum over several 
frequency bins. For example, the detection probability is higher than 90 % for modes 
spanning at least eight frequency bins with a signal-to-noise ratio greater than three. 
Therefore, we have searched the contemporaneous, very low-frequency {v < 500 /iHz) 
BiSOfJ^l. GOLF and SOI/MDI data for evidence of prominent peaks with various 
widths. We searched for clusters containing 2, 3, 4 and 5 spikes. We allowed a cluster 
to be spread over twice the width of the mode and took the widths of the modes to 
be the number of frequency bins covered if the mode lifetimes were 1 month, 2 months 
and 4 months, respectively (i.e., 32 bins, 16 bins and 8 bins, respectively). However, no 
statistically significant peaks were found within 4 ^iHz of the model g-mode frequencies. 



5.3 Innovative methods using a single set of data 
5.3.1 Time- distance analysis 



Duvalll |2004h devised a method for detecting g modes using time-distance helioseis- 



mology of p mode signals. The idea is to detect the small flows induced by the g modes 
on the travel times of the p modes. He used deep-focusing rays to create travel-time 
maps for detecting perturbations at about 100 Mm below the solar surface, and then 
applied regular spherical harmonics masks on these maps to obtain power spectra of 
these travel times. While this method has not yielded any potential detections, plenty 
of work is need to develop the technique. 



Birmingham Solar Oscillation Network, ICh aplin ct alj Jl996f) 



The quest for the solar g modes 



59 



5.3.2 Cross spectra of GOLF data 



Very recently, iGrec et"aL used about 13 years of GOLF data to search for g 



modes using cross spectra. The idea is to use about 22 sub-series, each of 7-month's 
duration, to create cross spectra that are calculated by pair-by-pair. In total, this 
provides more than 253 cross spectra (i.e., 22 x 23/2). The averaged cross spectrum is 
then normalised to remove the low frequency dependence, and then an autocorrelation 
of the cross spectrum is calculated. The autocorrelation of the cross spectrum shows 
a maximum at a value that is tentatively identified as being the g-mode rotational 
splitting. In addition, they also tried to detect in the 253 cross spectra the signature 
of peaks having the same phase. They compared their result with computed g-mode 
frequencies of a solar model, but found no convincing detections. 



5.4 Use of contemporaneous data 
5.4.1 A signal at 220.7 [iHz? 

At the beginning of the SOHO mission, there were several reports of the detection of 
a sign al close to 220.7 AtHz, that might possibly be identified as t he 1=2, n = —3 g 



mode dGabriel et al.ll 19991 ; Finsterle fc FrohlichlboOlh . Very recentlv jjimenez fc Garcia 
(2009]) and Garcia et al. ([2008a) reported again on the presence of this signal in data 



from the different intensity instruments comprising the VIRGO package on SOHO. 
There is a trace of the signal in GOLF (but not with the same quality), in GONG (but 
at the level of the noise) and nothing in SOI/MDI. However, it seems that the signal 
has no obvious instrumental or SOHO-related origin. It is possible that the origin of 
the signal is related to the noise created in intensity by convection close to a g-mode 
frequency, or in fact could be a component of a g-mode. 

5.4-2 Coincidence search 

Statistical tests have been used to search contemporaneous BiSON, GOLF and SOI/MDI 
data for low-frequency p modes, g modes and mixed modes using a frequentist ap- 
proach. To take advantage of the availability of contemporaneous data from different 
instruments frequency-amplitude spectra were searched for prominent spikes or pat- 
terns of spikes that were positioned in the same frequency bin or bins in any two 
of the three spectra. Details of the statistical tests can be found in iBroomhall et al] 
( 20071 ). Amplitude threshold levels were calculated for a 1 % probability of detecting a 



prominent feature by chance at least once in a frequency range of 100 /iHz. 

Figure [18] shows visually, in the form of an echelle diagram, the locations in fre- 
quency where prominent spikes or patterns of spikes were observed. Also plotted on 
Figure [TS] are the mode frequencies predicted by solar models. Detections were only 
considered as possible mode candidates if they were positioned within 1 /iHz of a pre- 
dicted mode frequency and such candidates are highlighted by a green square in Figure 

In the p-mode frequency ra nge all of the cand i dates correspond to previously 



claimed detections of p modes (See Toutain et al.l 19981; iBertello et ai]|2000al ; iGarcfa et al] 
200 lbl ; IChaplin et aflbooj IGarcfa et al]|2004 IBroomhall et al.ll2007I L Some of the de- 



tections in the g-mode range lie close to the predicted frequencies. However, it should 
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Fig. 18 An Echcllc plot, modulo 145 /iHz, marking locations in frequency of occurrences 
uncovered by the test searches. Locations in frequency where spikes, or patterns of spikes, 
were found in the same frequency bin, or bins, in BiSON and GOLF spectra at levels sufficient 
to record p < 1 % are marked by the black symbols in the middle of each row. A different 
symbol has been used for each test (see figure legend) . We have also recorded prominent spikes 
or patterns of spikes found by comparing either the BiSON and SOI/MD1 spectra (red symbols 
at top of each row) or the GOLF and SOI/MDI spectra (blue symbols at bottom of each row). 
The orange ve rtical lines mark locations of the frequencies of p modes predicted by the Saclay 
seismic model ((Turck-Chiczc ct al. 2 00lh . The vertic al purple lines mark locations of the m = 
frequencies of g modes predicted by the Ml model llProvost et al.|[2000l) . 



be noted that, since the number of g modes that lie in the frequency range searched is 
large, the chances of mistaking a false detection for a g mode candidate is large. There- 
fore, to improve confidence in the detections we required that each mode candidate 
passed more than one of the statistical tests, thus significantly reducing the number of 
false detections. This extra requirement discounts all of the detections in the g-mode 
range and some of the detections in the p-mode range. 

It is apparent from Figure [18] that the analysis has uncovered several occurrences 
of p < 1% that lie well away from the predicted mode frequencies. The number of 
these detections exceeds the number expected given the range in frequencies searched 
and the threshold probability. This possibly indicates that the statistical distribution 
in the noise is slightly different from the Gaussian distribution, which was assumed 
when deriving the statistical tests. Alternatively it could be due to the misleading 
significance levels assigned when adopting a frequentist approach (see Section 4.4.2). 
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6 Discussion and Conclusion 

Figure [15] shows the comparison of the best observational measurements with the- 
oretical amplitudes of the g modes. The theoretical g-mode amplitudes range from 
10 -2 mms -1 to a few mms -1 . The most optimistic theoretical mode amplitudes could 
be collectively detectable, after the collection of more than 10 years of observations of 
the solar radial velocity. Unfortunately, the detection of individual peaks is very far 
from being feasible since the only possibility would be to have the most optimistic 
amplitude larger by 50%, or twice as much in power (jBelkacem et alJl2009h . Since the 



noise reduction scales with the observing time T like log(T)/vT, another 80 years or 
so of data would be required to reduce the actual detection limit by a factor two. The 
limit does not scale like l/VT because the probability limit is kept constant in a given 
detection window w hile the number of frequency bins in that window increases like T 



( Appourchauxll 1 99a ) . At the time of writing, there is indeed a consensus amongst the 



authors of this review that there is currently no undisputed detection of solar g modes. 

Despite this state of affairs, the search is not over yet. What can be done to turn 
our fortunes around? First of all, some of the detection techniques presented in this 
review paper that rely on the theoretical knowledge presented in Sections 2 and 3. 
Several are yet to exploit fully the theoretical information: this calls for the application 
of full Bayesian statistical inference on the data (yet to be done) . 

The observation of signals from the the solar limb will hopefully offer another way 
of improving the potential visibility of the g modes. These observations will be carried 
out by PICARD, which is due to be launched in 2009. 

Because the solar atmosphere is relatively noisy, reductions of noise levels can 
potentially be achieved by making multiple observations that are sensitive to pertur- 
bations at different heights in the atmosphere. Here, one would seek to rely on changes 
in the coherence of the noise with height. The GOLF-NG instrument recently put in 
operation uses this princ iple for measuring solar radial velocities in 8 d ifferent heights 
of the solar atmosphere ( Turck-Chieze et al]|2006l ; ISalabert et alJl200aV 



Development of methods based on time-distance helioseismology offer promises not 
only for detecting directly the g modes, but also for probing the solar core by making 
observations from two widely separated point of view. The technique is called stereo- 
scopic seismology, and may provide the structure and the dyna mics of the core with 
a completely diffe rent approach with the Solar Orbiter mission ( Kosovichev fc Duvalll 
l200ilGizonll2006f ). 

Finally, we could try to detect the perturbation of the gravitational potential cre- 
ated by the g modes. For instance, the LISAF^I mission may provide a lower detection 



limit than the current classical helioseismic techniques. IPolnarev et al. I |2009h deduced 



a lower detection limit for LISA that is about a factor five lower in amplitude than 
the actual GOLF limit for modes of m — ±2 (see Figure I15|l . This limit would be 
low enough to d etect g modes if they had amplitudes as high as those predicted by 
iBelkacem et all (|2009h . The variations in t he potential can be also detecte d by using 



las er ranging, e.g the A STROD mission ( Nill2007l ; [Appourchaux et al.ll2009r i. As shown 



bv lBurston et all ( 20081 ). the capability of ASTROD would be sufficient to allow the 



detection of 



iKumar et aiT(|l996l) 



g mode s even if their amplitudes were to be as low as those predicted by 



The ESA/NASA gravitational wave laser interferometric space antenna 
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We conclude with an optimistic comment. It is superficially unfortunate that we 
may not yet have detected g modes in the Sun. On the brighter side, on which we 
always prefer to be, it shows that we have before us a greater challenge which will 
yield greater satisfaction when we overcome it. We are all now much more prepared to 
continue the search. 
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A Nonlinear amplitude limitation 

A.l Illustrative case of one-dimensional acoustic wave 

Consider a toy acoustic mode in a uniform medium in which Eulcrian perturbations, denoted 
by a prime, to pressure, p, density, p, and temperature, T, are related by the perturbed equation 
of state at first order: 

p' = c 2 p' + ST' . (45) 

Acoustic modes confined between two values, and 1, of a spatial variable x are slowly damped 
by viscosity and excited by a thermal process which we model crudely, and approximately, by 

8T' 

■w = -*> (46) 

where u is velocity. The following equation then governs the wave motion: 

5V a ay a f ay 

— o c 5 I *P "r v T I 

dt 2 dx 2 dt \ dx 2 J dx \ dx J 

where c 2 , p, A = S/3/cp and the kinematic viscosity v are all constants. We have also introduced 
the device of replacing d/dx by c~ 1 d/dt in the thermal term, which is now Xdp' /dt, to ensure 
that it excites both forward and backward propagating components of the mode. A more 
physically realistic model not requiring that artificial device for the excitation could easily have 
been adopted, but the governing equation would have been unnecessarily more complicated. 
The Eulerian density perturbation may be represented as a sum of normal modes: 



d 2 p' 2 d 2 p' d / , d 2 p'\ 8 ( du\ 



A n e-^ + A* n e i ^ t )cosk n x (48) 



with n = corresponding to the parent mode and n = 1, 2 to the daughter modes, where the 
complex amplitude A n (t) is slowly varying and the asterisk denotes complex conjugate. The 
wavenumber k n is determined by the boundary conditions, and the real frequency cj n satisfies 
uj n = ck„, so that the amplitude of a linearized unforced mode grows exponentially with time 
at rate tr„ = (\—uk^ l )/2, whose magnitude we consider to be small compared with ui n . It is now 
straightforward to show that when the weak nonlinear interactions between the parent and the 
two low-amplitude daughters are taken into account as a small perturbation, by substituting 
the linear eigenfunctions of Eq. H48H into Eq. H47I I. keeping only second-order terms, and 
integrating over x, the resulting equations for the amplitudes are given approximately by 



(49) 
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and 



dt 



■ 01,2-Ai 2 



-iAujt 



(50) 



in which I n = {p' n p'n) ls a measure of the inertia of mode, the angular brackets denoting 
integration over x, and A = [cos Akx]^ /AAk is a mode coupling constant which is an integral 
of terms trilinear in eigenfunctions; also -ui 2 are the intrinsic damping rates of the daughters, 
and are presumed to be positive, Au = u)q — w\ — u>2 is the frequency mismatch, and Ak = 
ko — ki — &2 dep ends on the ch ose n boun dary conditions. These equations are very similar to 
those derived by Dziembowski (sec 1982) for stellar g modes. 

Equation 1501 1 describes the pulchraf^l dynamics: if it is accepted that the daughters 
dissipate much more rapidly than the parent can grow, then the variation of Aq may at first 
be neglected; the equations are linear in j4„jq, and admit solutions A n = Q n exp(— Ldurt/2) 
with Q n oc exp(st) in which s > if \Aq\ 2 > A 2 , where 



Ill2U 1 U>2 



A 2 L 



1 + 



1x2 — <Ji 
02 + cr 1 



(Au,y 



(51) 



As we discuss later, A c estimates the limiting amplitude of the parent. Note that for a given 
coupling constant A and close resonance (Alo very small), the amplitude above which the 
parent excites her daughters is approximately proportional to the (harmonic) mean of their 
damping rates, and inversely proportional to the coupling constant. It is independent of the 
intrinsic growth rate of the parent (of course) , provided that the growth rate is small compared 
with the damping rates -cr n . A state of steady amplitudes exists when |Ao| = A c . In that case, 
other things being equal, the amplitude of the main mode would decrease if the viscosity were 
to be reduced, a result which at first sight might seem counterintuitive. 

To proceed further we need to apply these results to solar g modes, which requires consid- 
ering the g-mode frequency spectrum in spherical symmetry. Not surprisingly, for a nontrivial 
coupling constant one needs mi + 7712 = mo, and also 1% — 1\ = Zq or lo — 2, where m and I are 
azimuthal order and degree. These selection rules come f rom the integral of s pherical harmonics 
that appears as a factor of the coupling coefficient (sec Dziembowski 1982, for details) 



Yi Y 2 sin 



(52) 



where Yo,i,2 denotes the spherical harmonic associated to the modes (0,1,2), and (9,4>) are 
spherical polar angles. If the coupling constant is not to be small, one needs the orders m, 
n,2 not to be very different, so that the product of the daughter eigenfunctions varies on a 
scale similar to that of their parent. These conditions, together with the resonance condition 
\Aw\ « (^0; are most likely to be satisfied when l± and I2 are large, requiring I2 — li and 



consequently LU2 



jojoi because when £ is large the frequency distribution is dense, as 



can be seen from the approximate asymptotic eigenfrequency equation derivable from Eq. (2}: 



1/2 



d lnr : 



(53) 



the integral being between two radii at which N = u>. From that relation it follows that the 
frequency difference 8u between two modes of consecutive order satisfies Sui ~ fiojl" 1 with 
fi of order unity. Also the damping rates of the daughters satisfy —02 — —01, which, for 
convenience, we here denote simply by r). 



A. 2 Effect of frequency mismatches 

The scenario proposed by Dziembowski (1983) is that an 'equilibrium' state of steady oscil- 
lations is reached (ignoring temporal variation of the background state of the Sun) in which 



21 derived from Appius Claudius Pulcher, Roman politician of the 1st century BC who had 
several daughters and grand-daughters all named Claudia. 
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twin daughters of high degree extract energy nonlinearly from the parent and subsequently 
dissipate it at the same rate that the parent extracts en ergy linearly from the background 
state. Following a similar study by [Wersingcr ct al. (1980), Dziembowski also demanded that 
the equilibrium be stable, which requires Am > 1r\. Then, not knowing the structure of the Sun 
precisely enough to calculate the resonances, he considered instead the probable amplitude of 
the parent, making certain assumptions about the probability distribution of frequencies: the 
probability of being within Aw of resonance is proportional to Auj/Suj, and since 5u) ~ l~ 1 one 
might expect to find amplitude limitation to be most likely amongst daughters with greatest 
I. But also one must recognise that damping rates increase with I, according to r\ ~ rjol 2 (for 
a damping dominated by radiative losses), rendering the energy-extraction process at high I 
less effective and thereby permitting more freedom for the parent to grow. The outcome is 
a compromise. Dziembowski estimated the probability P that the parent has its amplitude 
limited to A to be 



the summation being over all pairs i of twins (having coupling constants Ai). Of course, the 
formula for the coupling constant, which depends on the precise way in which one chooses to 
define A, is much more complicated than that for the simple acoustic interaction discussed at 
the beginning of the section. 
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